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Abstract 



The crucial role played by pressure in general relativity is explored in the mathematically 
simple context of a static spherically symmetricc geometry. By keeping all pressure terms, 
the standard formalisms of rotation curve and gravitational lensing observations are extended 
to a first post-Newtonian order. It turns out that both post-Newtonian formalisms encode 
the gravitational field differently. Therefore, combined observations of rotation curves and 
gravitational lensing of the same galaxy can in principle be used to infer both the density 
and pressure profile of the galactic fluid, whereas the currently employed quasi-Newtonian 
formalisms only allow us to deduce the density profile. 

If a suitable decomposition of the galaxy model is used to separate the dark matter from 
the galactic fluid, the newly introduced post-Newtonian formalism might allow us to make 
inferences about the equation of state of dark matter. While the Cold Dark Matter paradigm 
is currently favoured in the astrophysics and cosmology communities, the formalism presented 
herein offers an unprecedented way of being able to directly observe the equation of state, and 
therefore either confirm the CDM paradigm or gain new insight into the nature of galactic 
dark matter. 

In a logically distinct analysis, I investigate the effects of negative pressure in compact 
objects, motivated by the recently introduced gravastar model. I find that gravastar like 
objects — which have an equation of state that exhibits negative pressure at the core of the 
object — can in principle mimic the external gravitational field of a black hole. Unlike a black 
hole, however, gravastars neither exhibit a pathological curvature singularity at the origin nor 
do they posess an event horizon. Instead they are mathematically well defined everywhere. 
Finally, another exotic option is considered as a mathematical alternative to black holes: The 
anti- gravastar, which is characterized by a core that has a negative mass-energy density. 
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Chapter 1 

Introduction 



Almost a century after Einstein's publications developing his general theory of relativity, it 
is astounding that the mathematics of the simplest case of static spherical symmetry still is 
not completely exhausted. In this thesis, I am going to discuss various aspects of general 
relativity that are of importance to today's research and that should be considered for the 
interpretation of current observations. 

Before I elaborate on these issues, I introduce the necessary mathematical framework of 
general relativity in §2. This includes mention of the basic ideas of general relativity, i.e. the 
equivalence principle, the equations of motion in general relativity, Einstein's field equations, 
and stellar structure in general relativity. I will also introduce curvature coordinates, which 
are mostly used throughout this work, and discuss the individual parts of the stress-energy 
tensor, as they are important for the following arguments. 

The first topic that I want to discuss in a general relativistic context is the gravitational 
field of galaxies. Today, the assumption that every galaxy has a dark matter halo is common- 
place, and usually justified by the observed rotation curves. However, the nature of this dark 
matter is completely unknown, and the only property that has been pinned down by observa- 
tions so far is the density of dark matter in galaxies. This density has been extracted from two 
kinds of observations: (i) galactic dynamics in which the gravitational field is inferred from 
the relative speed of involved particles, stars, gas, etc. and (ii) gravitational lensing where the 
deflection of light is used to deduce the gravitational field of the host galaxy. Since the grav- 
itational field of a galaxy is weak, both methods are usually analysed in a (quasi-)Newtonian 
approximation . 

I will show how to interpret both kinds of observations in a general relativistic weak field 
approximation that extends the usual quasi-Newtonian picture. As a consequence, combined 
measurements of galaxy dynamics and gravitational lensing can in principle yield additional 
information beyond the density of the galactic fluid, which consists of the visible parts of a 
galaxy and the dark matter. Combining both methods also allows one to deduce the pressure 
of the galactic fluid. If the accuracy of these measurements is sufficient, it may be possible 
to infer both the density and the pressure of galactic dark matter, and therefore gain infor- 
mation about the dark matter's equation of state. This could turn out to be unprecedented 
observational evidence for Cold Dark Matter (CDM), or even an indication of a new form 
of astrophysical field, such as e.g. Scalar Field Dark Matter (SFDM). Independent of what 
the actual result will be, in §3 I present a new formalism that is in principle compatible with 
current data reduction techniques and yields both the density and pressure of a galactic fluid 
from combined measurements of galactic dynamics (§3.2) and gravitational lensing (§3.3). 

The second aspect of spherically symmetric general relativity which I want to discuss in 
this thesis is a mathematical alternative to black holes. Just like dark matter, black holes 
are also widely accepted in the astrophysics community. Black holes are the most dense and 
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compact clumps of matter that can possibly exist. However, from a general relativistic point 
of view, black holes exhibit a pathological property: The mathematical solution of a black hole 
suggests that there is a curvature singularity at the center of the black hole, which basically 
translates to a failure of general relativity to describe this singular point appropriately. This 
singularity is only acceptable because it is hidden behind the event horizon, which is like a 
point of no return from which even light cannot escape. Therefore, one has no possible means 
(even in principle) of measuring anything behind that horizon. 

In §4, 1 will discuss these problems in more detail along with two mathematical alternatives 
to black holes that exhibit the same external gravitational field but do not show any of the 
previously mentioned pathological behaviour. The first alternative is the gravastar (§4.2), 
a concept that has been previously explored by Mazur and Mottola [58], or similarly by 
Laughlin et al. [19, 20]. I will extend their layered model to a continuous model and by that 
show that a basic property of a gravastar is the existence of anisotropic pressures, similar to 
a shear stress in an anisotropic crystal. The second alternative is the anti- gravastar model 
which was developed by Matt Visser and myself. For the anti-gravastar (§4.3), anisotropic 
pressure is not necessary, but one needs to accept a negative density at the core instead. Both 
alternatives are of mathematical nature only, and as of January 2006 there is no physical or 
observational justification to prefer them over the usual black hole picture. The discussion 
in this thesis merely illustrates that there are alternatives to black holes that do not exhibit 
pathological behaviour of the aforementioned kind at all. 

After the detailed dicussions in chapters §3 and §4, I will conclude in chapter §5 and point 
out the most important features of the new findings in this thesis. The following appendices 
contain several articles arising from the material in this thesis that have been published or 
submitted for publication. 



Chapter 2 

Mathematical framework of general 
relativity 

In this chapter I give ah definitions for the basic concepts and symbols I will use in later 
chapters. 

Geometric units are used by default unless otherwise stated. Thus, the speed of light 
c = 1, and Newton's constant of gravitation Gat = 1. Hence, all of mass, length, time, energy 
and momentum will have units of length [L\ and density and pressure will have units 
Multiplication by appropriate combinations of c and Gn yields Sl-units. 

The Einstein summation convention of omitting the sum symbol whenever a pair of 
contravariant (index up) and covariant indices (down) appears in one term is used throughout 
the whole thesis: 

3 3 

^M^'^^E^M^^ x.y^'^^x^y". (2.1) 

/i=0 a=l 

Greek indices indicate four spacetime dimensions . . = 0, 1, 2, 3) and latin indices indi- 

cate three spatial dimensions (a, 6, ... = 1, 2, 3). 

In places where the dimensionality of the index is unambiguous, I sometimes use the 
"bullet notation" to reduce the number of lettered indices and therefore increase the legibility. 
Bullets ("•" and "o") indicate indices that are used for contraction: 

A,B' = A^B^'- A.O B^'* = A^p B^''' . (2.2) 
2.1 Equivalence principle 

The theory of general relativity is based on the principle of universal free fall, also called the 
(weak) equivalence principle. From the observation that the inertial mass and the gravita- 
tional mass are identical to remarkable precision, one can conclude that all objects "fall" in 
the same manner, i.e. the actual trajectory of a body under the influence of a gravitational 
force is independent of its mass. This is not true for other forces like e.g. the electromagnetic 
force. 

Combined with the special theory of relativity, Einstein formulated the strong equiva- 
lence principle, also called the "Einstein equivalence principle", which states that gravity is 
represented by the Christoffel connexion on a Lorentzian manifold with an associated metric 
tensor g^^y, and that: 
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• the universality of free fall is given by the geodesies (2.5) of the Christoffel connexion 
(2.9), and 

• the rules of special relativity are recovered in a local rest frame, where the metric tensor 
takes Minkowskian form (2.12). 

This approach guarantees that the equivalence of masses is imposed at a very basic level. 
Since the equations of motion are given by the Christoffel connexion and hence, the metric 
tensor, Einstein still needed to find field equations (2.33) that relate the metric to the mass- 
energy which acts as the source of the gravitational field. 

I will now summarize these fundamental aspects of general relativity to set the framework 
for the following chapters. 



2.2 Geodesies and affine parameters 

Geodesies are defined as the "straightest possible" curves in a manifold. Straight means that 
tangent vectors to a geodesic remain parallel to each other along the curve. Parallelism in 
a manifold is given by the notion of an affine connexion T^^^. Please refer to the standard 
literature about differential geometry and non-Euclidean spaces for further details, e.g. [35, 
60, 89, 92]. 

The covariant derivative V/j with respect to the coordinate is the generalization of 
taking derivatives in curved spaces. The departure from flat space is represented by the affine 
connexion P'^^^. Applied to a vector A", the covariant derivative is given by 

VpA^^dfsA'^ + T^pA', (2.3) 

where 9^ denotes the partial derivative with respect to the coordinate . Without deriva- 
tion, I note that a curve with an arbitrary parametrization x is a geodesic iff its tangent 
vectors, 

dX^ 

t'^^, (2.4) 

satisfy the differential equation 

^+P- ^^-f(x)^ (25) 

or equivalently, 

^^v^^'^ = /(x)^^ (2.6) 

where /(x) is an arbitrary function of the parameter x- If X is chosen in such a manner that 
fix) vanishes, x is called an affine parameter. In that case it is possible to write (2.6) in 
index free notation, 

t-Vt = 0, (2.7) 

to illustrate that the tangent vectors t along a geodesic curve have constant length. This 
is exactly what an affine parameter does, it parametrizes the curve in such a way that the 
"parameter speed" along that curve is constant. For later reference, I note that in the usual 
component notation, the geodesic equations for an affine parameter x, (2-7), take the form 

d^Xt" dX" dXf^ _ 
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2.2.1 Metric spaces 

The metric connexion without torsion is the usual connexion that is used in general relativity. 
In this case, the affine connexion is given by the Christoffel symbols of the first kind, 



^tiap = {a*; a/3} = ^ {da 5m/3 + df3 gfj,a - d^Qap) 



(2.9) 



or the Christoffel symbols of the second kind. 



0/3 



a(3 



^ 9^' {da g,f3 + 0/3 g,a - d,gap) 



(2.10) 



where the contravariant metric tensor g'^'^ was used to raise the index To lower an index 
one uses the covariant metric g^^: 



A. 



g 



(2.11) 



In flat space and local orthonormal frames, the Minkowski metric is used to raise and lower 
indices: 

??M>^ = r?^" = diag[-l, 1,1,1]. (2.12) 

For more details on dual spaces and index gymnastics see the introductory literature to 
differential geometry and general relativity, e.g. [35, 60, 89, 92]. 

In metric spaces the geodesic notion of "straightest possible" curve coincides with the 
shortest possible curve in Riemannian geometries (with positive definite metric tensor) and 
with an "extremal distance" in Lorentzian geometries (where the metric tensor has the sig- 
nature — h + +, i.e. that of the Minkowski metric ^^v). 

Let gp,y be the metric of the 3- or 4- (or even n-) dimensional Riemannian or Lorentzian 
space, and x be a parameter of an arbitrary curve. Then I define arclength s of that curve 
analogously to Euclidean arclength in flat space [89], 



ds 



dx = 




dXM AX^ ^ 
gfiu — — — dx • 



dx dx 



(2.13) 



Applying the Euler-Lagrange equations of variational calculus yields a differential equation 
that characterizes the path which has the shortest arclength s between the fixed endpoints 
of the curve along which the integral is evaluated: 



dx 



d{dX-/dx) 



d 



9fiu dx dx 



0, 



(2.14) 



which is 



d 

dx 



■ga/s 



dX" 



^ dxt^ AX'' ^ dx 
9iiv dx dx 



1 



dga^ dX<^ dX^ 

^^I9ixv dx dx 



(2.15) 



If the curve is not a null curve in a Lorentzian space (for which s = everywhere), this can 
be reparametrized in terms of the arclength s. Hence, I multiply by 
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which simpHfies (2.15) tremendously to 



9al3 



ds 



1 dgp^ dX^ dX^ 

2 dX^ "dT "dT 



Expanding the first term and rearranging the derivatives of the metric gai3 gives 



gal3 



d^ , idg^f, 1 

ds2 



+ 



dXf^ dX^ 



dX^ 2 ds ds 



(2.17) 



(2.18) 



which is easily identified as the geodesic equation with the ChristofFel symbols of the first 
kind (2.9), 

d^X^ ^ dX^ dX^ 

+ -T- -J- = ' (2-19) 



or equivalently, using the Christoffel symbols of the second kind 



d^ 

ds2 



+ r' 



dXf^ dXT 
ds ds 



0. 



(2.20) 



Hence, a geodesic describes the shortest possible path in a Riemannian space where arclength 
s is defined as in (2.13). In a Lorentzian space, this calculation requires ds 7^ 0, hence it does 
not apply to null curves. For all timelike and spacelike curves, the geodesic is equivalent to 
the curve with a locally extremal arclength s. 

This calculation also shows that for metric spaces, arclength s is an affine parameter, 
since (2.8) is always satisfied for positive definite metrics or spacelike and timelike geodesies 
in Lorentzian spaces. 



2.2.2 AfRne parameters for null curves 

The affine parameter for lightlike geodesies, also called null curves, cannot be easily specified 
in general, since the motion of massless particles at the speed of light is characterized by the 
vanishing of the invariant interval, 

ds^ = g^^ dX^' dX" = . (2.21) 

Hence, the previous derivation of affine parameters for time- and spacelike curves is invalid 
for null curves, due to a division by zero. Instead the affine parameter has to be determined 
for each metric tensor g^y. I will only consider static spacetimes, which are the only kind of 
spacetimes used in this thesis. 

If a spacetime is stationary, the metric is independent of the time coordinate, dt g^u = 
and therefore, by the invariance of infinitesimal coordinate translations in the t-direction, 
K^^ = 5^^ is a timelike Killing vector [60, §25.2]. From the Killing equation follows that the 
covariant derivative of the Killing vector is completely antisymmetric in its two indices and 
hence, 

fV,{fKo) = {fV.t°)Ko + ff{V,Ko) = {t'V,t°)Ko , (2.22) 

where t'^ = dX^/dx is tangent vector to a null geodesic. But by (2.8), this must vanish if x 
is an affine parameter. Consequently, the quantity 

. n dX' dt f dX^\ , , 

E = fK, = g,o fK° = g,t —— = — [ gu + gu -r- ] (2.23) 



is conserved along all geodesies with an affine parameter - it does not matter whether the 
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geodesic is null or not. However, since the affine parameter for null geodesies is still lacking 
an operational definition, I rearrange (2.23) into a definition for the affine parameter x of 
null geodesies, 

(dX* \ 
gtt + gu^j , (2.24) 

where the constant —E was absorbed into the affine parameter x, as the overall normalization 
of an affine parameter is irrelevant. In the more restricted case of a static spacetime (with 
gti = 0), the affine parameter of null curves takes an even simpler form: 

dx = -dtgtt. (2.25) 

As a note on the side, this can also be expressed using the three-dimensional proper distance, 



dx = V^t ^9^j dXi , (2.26) 
which can be derived with the help of the invariant interval (2.21) for static spacetimes: 

gttdt^ + gijdX'dP = Q. (2.27) 



2.3 Curvature coordinates 



In this work, I will concentrate on static, spherically symmetric fluid spheres. It is a standard 
result, see e.g. [60, §23.2] or also [92, §8.1], that all geometries of this category are represented 
by a metric of the form: 

ds^ = g,,u dX^ dX"" = -A{r) dt^ + B{r) dr^ + d^J^ , (2.28) 

where g^u are the metric components and X^^ is the set of coordinates. The arbitrary functions 
A{r) and B{r) depend on r only^, t is the time coordinate, r is the space coordinate in the 
radial direction of the sphere and 

dO^ = d6'2 + sin^ (9 d(^2 (2.29) 

is the metric of the two-dimensional unit sphere with the two spherical polar coordinates 9 and 
Coordinates of this type are called "curvature coordinates" or sometimes "Schwarzschild 
coordinates" . This choice of r coordinate has the advantage of a physically clear meaning of 
r = (proper circumference) /27r, where "proper circumference" refers to the integrated proper 
distance interval, ds, on a closed circle with its center at the coordinate origin: 

(proper circumference) = / ds = r dip = 2'nr. (2.30) 

Jip=0 Jo 



^Weinberg [92, §8.1] calls metrics of this type "static and isotropic", where static refers to independence 
of the time coordinate (as usual) and isotropic is the property that the metric is invariant under rotation 
- which is identical to the term spherically symmetric in this thesis. Weinberg's slightly unusual use of the 
term isotropic is not to be confused with isotropic fluids that will be mentioned later on, nor with the distinct 
concept of an isotropic coordinate system which is also commonly used in general relativity, and later in this 
work. 
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Here dt = 0, dr = and, due to spherical symmetry one can choose 6 = 7r/2 [60, §23.3]. 
While this seems to be trivial at first, note that the proper distance in radial direction is 

/•r=r2 rr2 

(proper radius) = / ds = \J B{r) dr / r2 — ri . (2.31) 

Jr=T\ J ri 

Depending on the problem under investigation, it is generally advisable to replace the 
functions A{r) and B{r) by functions which are more descriptive in the context of that 
problem. I choose to use A{r) = exp[2<I>(r)] and B{r) = 1/(1 — 2m{r)/r) throughout this 
work unless otherwise stated: 

= -e2*W dt^ + l!i^ + ^2^Q2_ (2.32) 

1 — 2m(r) /r 

To understand the physical meaning of the two metric functions, $(r) and m{r), it is necessary 
to invoke the Einstein field equations. 



2.4 Einstein equations 

The Einstein field equations are derived in basically every standard textbook concerning 
general relativity, for example [35, 60, 92]: 

Rfj-u - ^9iiu R = 8tt T^j, (2.33) 

where Rfj,u is the Ricci curvature tensor, R = g'^'^ R^i, is the Ricci scalar and T^u is the 
stress-energy tensor. In Sl-units, Svr is replaced by SttGn/c^. The combination 

Guv = R/iu — ij^a^iv R (2.34) 

is called the Einstein tensor, which can be obtained through a tedious but straightforward 
calculation once the metric components g^j_i, are given. Nowadays, this monotonous task is 
usually executed by a computer algebra program such as Maple^ or Mathematica'^. 



Geometric part of the field equations. The left hand side of (2.33) represents the 
geometry of the spacetime and is given as a non-linear combination of the metric components 
g^u and their first and second derivatives. For completeness, I shall define the Riemann 
curvature tensor, 

= 9« r% - dp r^,„ + r^.„ r',^ - r'^.^ , (2.35) 

and the Ricci tensor which is given by the contraction over the first and third index of the 
Riemann tensor: 

R^u = R*^i,v ■ (2.36) 

As a purely geometric result of metric spaces, the Einstein tensor obeys the contracted Bianchi 
identity, 

G"^ = . (2.37) 



© by Maplesoft, a division of Waterloo Maple Inc. 
© by Wolfram Research, Inc. 
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Components of the Einstein tensor for given metric. Since the Einstein tensor is 
completely determined by the metric, I shall give the components of the Einstein tensor, 
as derived from (2.32), for further reference. The significance of second rank tensors with 
one index up and one index down is explained later in §2.5.2. Prime denotes the partial 
derivative with resepect to the r-coordinate: ' = d/dr. The follwing components of the 
Einstein tensors are the only non-zero components of the metric (2.32) which were derived 
using a computation by Maple and manipulations by hand: 



r' 



G\ = 



2 m' 



m , I 2m 
r <^ 1 



(m' r — m)(l + r $') 



+ 1 



2 m 

r 



r 



(2.38) 
(2.39) 

(2.40) 



2.5 Stress-energy tensor 

The righthand side of (2.33) is the source term in form of a stress-energy tensor that describes 
the matter or field which is creating the curvature in spacetime. The stress-energy tensor 
— also called more accurately the energy momentum tensor — is generally the sum of all 
different kinds of stress-energy present in the system under investigation, e.g. fluid stress- 
energy, electro-magnetic stress-energy, stress-energy arising from other fields, etc. : 

= + + Tfi^eids + • • • • (2-41) 

The stress-energy tensor of a swarm of n particles indexed by the particle number i is given 
by [92, 2.8.5a] 

n « 

VlrmiXn = E "^o,i / Vr Vf S\X>^ - Xf (r)) dr (2.42) 
1=1 

where 5^{-) denotes the four-dimensional delta- function and mo,j is the rest-mass of particle 
i with the four- velocity V"' = dX"/dr at the position given by the coordinates X-^. Fur- 
thermore T is the proper time that parametrises the trajectories of all particles. Instead of 
referencing the individual particles by their indices i, one can write the stress-energy tensor 
in terms of the local number density n(X^), the mass of the particle at X^, m{X^), and the 
local four-velocity y°(X^') of the particle at Xi" [35, 22.17]: 

TZi,^{X'^) = mnV'^V^. (2.43) 

At this point some amount of averaging takes place: the mass of the particle at X^ is actually 
the average mass of the particles in the immediate vicinity of X^, using the number density 
of the associated species as a weight for averaging. If one is only dealing with one species of 
particles in the fluid, m{X^) is of course a constant. 

From now on, I will drop the index "swarm" of the stress-energy tensor, since it is the 
only type that is considered in the rest of the thesis. However, I will use it in a slightly 
different form that corresponds to a continuous fluid. 

To see that, let us consider a small volume /NV at rest in the observer's orthonormal 
frame which is part of a constant-t three-surface with an associated timelike normal vector 
vp = (1;0;0;0). The product [35, 22.13] 

AP" = T^^ vp Ay = T"* Ay (2.44) 
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is identified with the four- momentum = nioV^ of the fluid in that smah volume [60, §5.3] 
so that we can attribute physical meaning to the components of the stress-energy tensor in 
the observer's orthonormal frame (denoted by the hats on the indices): 

r** = = p (mass-energy density) (2-45) 

T** = ~ '^^ (momentum density in direction i) . (2.46) 

Similarly, we can look at a three-volume AiAX^AX^ with a spacelike normal vector in the 
direction [35, 22.18]: 

AP" = T"^ AtAA^ (2.47) 

where IS. AS = AX^AX^ is the spacelike two-area normal to the direction. We can repeat 
this procedure for all three space directions and find: 

AP* AP* A.X^ 

T*-^ = .... = — : — = TT^ (energy flux in direction 7) (2.48) 

AA^At AV At V Bj jj \ J 

T^i = = = P (stress tensor) . (2.49) 

AA^ (area) 

Due to the equivalence of mass and energy in general relativity, one realises that momentum 
density and energy flux are actually the same physical quantity vr*. This has to be so since 
the stress-energy tensor is symmetric in its two indices, as can easily been seen from (2.42) 
and (2.43). The components of the 3-dimensional stress tensor f^^ are the i components of 
a force per unit area - also called a stress in classical mechanics - exerted across a surface 
with normal in direction j [35, 22.22]. 



The easiest example of a stress tensor is pressure in a fluid as seen by a comoving observer: 
the fluid is at rest relative to the observer and thus, the forces exerted is constant in all 
directions and the force per unit area is simply the pressure p [35, 22.25]: 

T^'^=p6'\ (2.50) 
where 6^^ = diag[l, 1,1] is the Kronecker delta with two contravariant indices. 



Again, some more averaging took place in going from a particle swarm to the notion 
of a continuous fluid. When we are talking about a force that is exerted across a surface, 
we need to realise that in a fluid, that force is created by momentum transfer between the 
interacting particles. So we either keep track of the exact motion of all particles, as is 
suggested by (2.42) and even (2.43), and the stress tensor drops out implicitly; or we average 
over the small fluctuating motion that creates the pressure (Brownian motion) and insert 
the stress tensor explicitly into the stress-energy tensor. Due to quantum fluctuations, there 
will always be some amount of random motion in the fluid. If one wants to speak about a 
fluid at rest - which is really a fluid in dynamical equilibrium - one has to average over the 
small fluctuations and introduce the pressure explicitly. If there is no interaction between 
the particles whatsoever, no pressure can be exerted within the "fluid". 



The notion of pressure is important in this context because it contributes to generating 
the gravitational field as I will soon show. This contribution arising from the pressure will 
play a vital role in the formalism which I will introduce in §3. 
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2.5.1 Perfect fluid 

An isotropic fluid, like the one that was considered in (2.50), is said to be perfect when heat 
conduction, viscosity or other transport and dissipative processes are negligible [35, §22.2]. 
Such a fluid is free of shear stress in the rest frame which implies that the stress tensor is 
diagonal and has 3 identical eigenvalues [60, §5.5]: 

oc 5'^ . (2.51) 

Hence, a perfect fluid in an orthonormal rest frame takes the form 

T^'^ = diag[p,p,p,p] = p6^^6\+p5^^5''/^ , (2.52) 

which can also be written in terms of the flat space Minkowski metric of the orthonormal 
frame, rj^'^ = diag[— 1, 1, 1, 1], and the four-velocity of the observer at rest = (1; 0; 0; 0): 

= {p + p)u''u'' +pr]f"' , (2.53) 

The natural generalisation of this perfect fluid stress-energy tensor is to permit arbitrarily 
moving observers with corresponding four-velocity u°' and to replace the locally flat spacetime 
metric of the orthonormal frame, tj^'^ , with the metric of the curved spacetime, g^'^ (we can 
now drop the hats, since we are no longer in orthonormal coordinates) [35, 22.39]: 

T^^ = {p + p) u^u" +pg^''' . (2.54) 

We can see that this generalisation is consistent with the Einstein equations (2.33) by noting 
that because of the general relativistic Euler and continuity equations, (2.54) satisfies the 
covariant conservation of stress-energy, 

V^r^'^ = 0, (2.55) 
which is necessary to fulfill the contracted Bianchi identity (2.37). 



2.5.2 Diagonal metric 

When working with a diagonal metric g^^^ one can benefit from using second rank tensors 
with one contra- and one covariant index. A short calculation shows that the components of 
a diagonal second rank tensor with one index up and one index down do not change when 
the coordinates are transformed into an orthonormal basis, i.e. 

A^^^ = A-p6\6^^^ = A-p. (2.56) 

An orthonormal basis (e^)" in a space with signature — h + + is characterized by the con- 
ditions [35, §7.8] 

go.(3{e&r{epf = ri^^ (2.57) 

where the indexing convention has to be understood in the following way: Indices without 
a hat label components in the curved space, the realm of gap. Indices with a hat indicate 
components of the flat space (with metric 77^^) that is tangent space to the curved space. The 
index a in (e^)" is a labeling index that enumerates the vectors in the basis while the index 
a is a coordinate index that refers to the basis-vector's a-coordinate in the curved space. 
If the metric gap is diagonal, finding an orthonormal basis is particularly easy since the 
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basis vectors of a.re already orthogonal. The only step left to do is to normalize the basis 
vectors at every point in the spacetime: 



(eaT = ; = W ^ K ■ (2.58) 

V 9aa y 9 

The transformation of a second rank tensor with one contra- and covariant index each is then 
given by 

A% = Ay {e% {e^r = A^ ^ 5/ ■ (2-59) 

If the tensor AP^^ is diagonal, the square root in (2.59) equals unity for every non-vanishing 
component and (2.56) is obtained. If the signature of the tensor A with both indices either 
up or down is assumed to be + + + + , then it follows from 

A'^=r?,^A-=r?-A^, (2.60) 

that the signature of the tensor A with one index up and one index down is — h + + . Hence, 
the stress-energy tensor (2.52), for example, takes the form 

= diag[-p,p,p,p] =r<|. (2.61) 

The minus sign in front of the density p does not indicate a negative density, but is a mere 
reflection of the tensor's signature — h ++. This slight awkwardness is the price for the 
complete absence of physically irrelevant coordinate artefacts in the components of the tensor. 



-p 
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2.5.3 Canonical forms of the stress-energy tensor 

The stress-energy tensor of the perfect fluid was one particularily simple example of how 
matter is represented in general relativity. In general, the stress-energy tensor has more than 
two independent components. Hawking and Ellis [38, §4.3] distinguish four possible canonical 
types of stress-energy in an orthonormal basis. The following classification is more or less a 
direct quotation from [38]. 



Type I is the general case where the stress-energy tensor has one timelike eigenvector that 
is unique (unless p = — pj). It can be expressed as 



p 














V\ 




























P3 



(2.62) 



The eigenvalue p represents the energy-density as measured by an observer at rest, i.e. the 
observer's four-velocity is u^y^^ = (1;0;0;0) in the orthonormal basis. The eigenvalues pi 
represent the principal pressures in the three spacelike directions. All observed fields with a 
non-vanishing rest mass are of type I, as are all zero rest mass fields with the exception of 
the special cases that are represented by stress-energy tensors of type II. 
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Type II is the special case where the stress-energy tensor has a double null eigenvector: 

(2.63) 



U + K V 

V V — K 

Pi 

p2 



The only observed occurrence of this form is for zero rest mass fields that represent radiation 
that is travelling only in the direction of the double null eigenvector. 



Type III &; Type IV are the cases where the stress-energy tensor has a triple null eigen- 
vector and no timelike or null eigenvector. Hawking and Ellis [38] are not aware of any 
physical occurence or relevance of types III and IV, and neither am I. 

In this thesis I will only focus on stress-energies of type I which represent the majority of 
relevant matter distributions. For simplicity, I will further restict the contents of this thesis 
to spherically symmetric cases. 

The nature of the involved pressure terms is irrelevant for the discussion presented. That 
is, it does not matter whether the pressure or stress of that stress-energy arises from the 
random motion of particles or from a field of some sort (see e.g. §3.1.5). 



2.5.4 Spherically symmetric non-perfect fluids (anisotropic fluids) 

In spherically symmetric coordinates, the rotation invariance of the sphere ensures that the 
physical properties of the 9- and (^-directions are identical. Mathematically this is established 
through the occurence of the unit sphere's metric (2.29) in the metric of the spacetime (2.32). 
Hence, the 69- and (/9(^-components of a spherically symmetric second rank tensor, that is 
derived from the metric and its derivatives, must be identical apart from coordinate artefacts. 
Thus, the Einstein tensor G"^ has identical 99- and (^93-components, as it is already evident 
from our example (2.38) to (2.40). 

Since the Einstein tensor is by (2.33) related to the stress-energy tensor, the latter one 
must also have identical 99- and c^yj-components to represent a matter distribution that 
sources a spherically symmetric gravitational field. Generally, any second rank tensor that 
decribes any spherically symmetric property must have identical 69- and 93(^-components. We 
can easily check this fact by noting that the components of such a tensor in an orthonormal 
frame must be invariant under rotation of space about the r-direction. Such a rotation is 
given by the rotation matrix 



n(7) 



1 










COS 7 
sin 7 






— sin 7 
cos 7 



(2.64) 



Rotating the tensor A'^g = diag[— a, b, c, d] by an angle 7 gives 



-a 

6 

0c cos^ 7 + (i sin^ 7 {c — d) sin 7 cos 7 

(c — fi) sin 7 cos 7 c sin^ 7 + c? cos^ 7 



(2.65) 



so that rotation invariance is only given when c = d. 

Hence, a stress-energy tensor that describes a spherically symmetric matter distribution 
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as measured by an observer at rest is given by 



-/> 

p,. 

pt 





(2.66) 



where pr and pt are the pressure or stress in the radial and transverse direction. 

Please note that "an observer at rest" is equivalent to saying that the observer's four- 
velocity has only components in the time direction, u"^^ = (1;0;0;0). Therefore, (2.66) is 
still a good approximation when the observer's three-velocity relative to the fluid's motion is 
small compared to the speed of light. 



The stress-energy tensors in this thesis are of the form (2.66), which represents the 
whole class of type I stress-energy of the Hawking and Ellis classification [38] subjected to 
spherical symmetry. This includes ordinary baryonic matter as well as some dark matter 
candidates and scalar fields as they appear later in §3.1.5. 

All static spherically symmetric fields and stress-energies of type I can be reinterpreted in 
the terminology of energy-density p and principal pressures pr and pt, and thus, are subject 
to statements made in this thesis. 



2.6 Stellar structure equations 

The structure of any static spherically symmetric system that obeys the rules of general rel- 
ativity is given by the "stellar" structure equations, see for example [60, §23.5]. The name 
stellar structure equations is commonly used since these equations are generally used to de- 
scribe compact objects like ordinary fusioning stars and neutron stars, etc. They do, however, 
also apply to larger systems like galaxies (in the spherically symmetric approximation). 

The structure is determined by five functions of the radial coordinate r: the metric 
functions $(r), m{r), the density p{r) and the radial and transverse pressures pr('^) andp((r). 
Therefore, to uniquely determine the structure, an appropriate model needs to supply five 
equations plus boundary conditions. The Einstein equations (2.33) provide three equations 
in static spherically symmetric coordinates: The tt-field equation yields 

/•r 

m{r) = 47r / p(f) df -|- m(0) (2.67) 
Jo 

and the rr-field equation gives 

= ^'(^^ = m{r)+4TTpr{r)r^ ^ 47rr p{r) + 3pr{r) 
y^')--^^') r^l-2m{r)/r] 3 l-2m{r)/r ^ ' 

where the "average density" p{r) = m{r) / {^irr^) was introduced. Please note that r is just 
a coordinate for the radial direction and not "proper radius"! The name "average density" 
is thus somewhat misleading, but very illustrative. 

From now on, in the interests of legibility, I discontinue indicating the explicit r-dependence 
of all relevant functions. The last non-zero Einstein equation could be used as a structure 
equation, but its complicated form is highly inconvenient. Indeed, one generally replaces the 
remaining field equation by the covariant conservation of stress-energy (2.55). To see that 
this is a valid way to proceed, let us count the equations at hand: The 10 components of the 
Einstein tensor by construction obey the four contracted Bianchi identities (2.37) and con- 
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sequently imply the covariant conservation of stress-energy (2.55). Thus, one can choose to 
disregard up to four Einstein equations by demanding that the Bianchi identities - or equiva- 
lently the covariant conservation of stress-energy - are also satisfied. The Einstein equations 
supply 10 equations, of which the six equations corresponding to the off-diagonal elements 
of the stress energy tensor are vacuous ("0 = 0") in the employed coordinate system with 
the diagonal metric (2.32). The two equations corresponding to the 6- and (/^-coordinates are 
identical due to spherical symmetry. The remaining three equations automatically satisfy the 
covariant conservation of stress-energy, 

p'^ = -ip + p^)^'+ ^^P'-P^\ (2.69) 

r 

and hence, we can choose to replace the 00-field equation with (2.69). 

The physical interpretation of p, pr and pt was already discussed in §2.5. The interpreta- 
tion of m(r) is given by (2.67): m{r) is the "total mass-energy inside the coordinate radius r" 
and this includes contributions from "rest mass-energy" , "internal energy" and "gravitational 
potential energy". For details, please refer to [60, Box 23.1]. 

To find the meaning g = <^', I point out that the gu component of an accelerating 
observer's proper reference frame is generally given by [60, eq. 13.71] 

gu = -(l + 2a.xJ) + o[|x^f] (2.70) 
= -{l + 2ajx^) + 0[\x^^] , (2.71) 

where is the physical acceleration felt by the observer. For comparison, the gu component 
of the metric (2.32) in the vicinity of tq is to first order in r 

g,^ = _e2i'(-o+5r) ^ _g2$(ro) ^ 3 «>'(ro) 6r] + O [r2] , (2.72) 

which leads, in a properly normalized reference frame with origin at r = rg, to 

g^iro + 5r) = - [1 + 2 ^'(ro) <5r] + O [r^] . (2.73) 

One can now easily identify 5 = <!>' = Or as the local gravitational acceleration felt by an 
observer in a proper reference frame, which is — for positive g — pointing towards the coor- 
dinate origin in the Newtonian picture and away from the origin in Einstein's interpretation. 
In analogy to Newton's gravity, we shall call <I>(r) the potential or even the gravitational 
potential. 

So far I have given you three of the necessary five structure equations and physical inter- 
pretations of all involved functions. The two remaining equations are given by the properties 
of the matter or field in our system. They usually occur in the form of equations of state 

Pr = Pr{p,r,n,s,...) (2.74) 
Pt = pt(.p,r,n,s,...) (2.75) 

but can they can generally also be replaced by other equations, see appendix §B for a more 
extensive discussion. I indicated that the equations of state generally depend on a variety of 
physical properties, like density p, number density n, entropy per particle s, or also chemical 
composition X, or even position r, etc. For each new function that adds complexity to the 
system, a further equation is necessary to determine the system uniquely. Within the scope of 
this thesis, I am only going to explore systems that are determined by the five aforementioned 
functions. 
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2.7 The Tolman-Oppenheimer-VolkofF (TOV) equation 

In most discussions, the derivative of the potential ^' is ehminated from the structure equa- 
tions by inserting (2.68) into (2.69). Thus, one obtains the anisotropic TOV equation: 

dpr {p + pr) {m + Airpr r^) 2{pt-pr) 

IT = rn — ^ — n + ■ (2-76) 

dr [1 — 2,m/r\ r 

In the case of isotropic pressures p = Pr = Pt this leads to the (standard) TOV equation, 
which can be written in any of the equivalent forms 

dp _ {p+p){m + 4Trpr^) _ Air r {p + p) {p + 3p) 
dr [1 — 2m/r] 3 1 — 2m/r 

To see the connection to Newton gravity, it is helpful to write the isotropic TOV in the 
following way: 

dp pm^^^pX / ^47rpr3\ / _2m\-i_ 



dr y P J \ 

This is the Newtonian equation of hydrostatic equilibrium multiplied by three relativistic 
corrections that are negligible in the Newtonian limit where |p| <C p and m <ti r. The inter- 
esting fact about this formula is that in Einstein's gravity it drops out automatically through 
the interaction between stress-energy and geometry while in Newton mechanics hydrostatic 
equilibrium is a seperate concept of fluid mechanics. Of course this only works because the 
appropriate fluid mechanics are already implemented in the way the stress-energy tensor was 
defined ~ in Einstein's theory of gravity, it is not possible to separate fluid mechanics and 
gravity! 



Chapter 3 

Galactic Dark Matter halos 



3.1 Gravitational field of a galaxy 
3.1.1 Types of galaxies 

A galaxy is a gravitationally bound system of lO^*' to lO"*^^ stars [8, §1.1]. According to 
Hubble's revised classification of 1936 [42], one distinguishes four basic types of galaxies, 
elliptical (E0-E7), spiral (Sa-Sd), spiral bar (SBa-SBd) and irregular (Irr), and the transition 
class of lenticular galaxies (SO). 




Figure 3.1: Hubble's classification of galaxies . 

Elliptical galaxies are smooth and featureless stellar systems that typically consist of old 
stars with low heavy metal abundances (Population II stars). Elliptical galaxies contain little 
or no interstellar gas or dust. In the classification, they are denoted by En where n is ten 
times the ellipticity of the galaxy, thus EO being spherical and E7 being highly elongated 
galaxies. [8, 61] 

^Figure taken from http://hubblesite.org, reproduced with permission of NASA and STScI. 

1 1 
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Lenticular galaxies are the class in transition between the elliptical and spiral galaxies. 
They are smooth featureless disks without gas, dust or bright young stars. They do obey 
the exponential surface brightness law of spiral galaxies and are labelled SO in Bubble's 
classification. [8] 

Spiral galaxies, denoted by Sa-Sd, contain a rotating disk of relatively young metal-rich 
stars (Population I), gas and dust. Within the disk are dense filaments of bright young stars, 
gas and dust that demark active star formation areas. These filaments take the shape of two 
or more spiral arms and are responsible for the name of this class. Apart from the disk, spiral 
galaxies also contain a central bulge and a spheroidal halo of Population II stars that form 
globular clusters. Along the sequence Sa-Sd, the relative luminosity of the bulge decreases, 
the fraction of gas increases, and the arms become more loosely wound. The Milky Way is 
classified as Sbc, an intermediate of Sb and Sc. [8, 61] 

Spiral bar galaxies are distinguished through the existence of a central bar from which 
the two main spiral arms begin. A bar is effectively an elongated central bulge. Apart 
from the added bar, spiral bar galaxies follow the sub-classification of normal spirals and are 
distinguished by the prefix SB. [61] 

Finally, the class of irregular galaxies (Irr) comprises all other galaxies that don't fit into 
the other classes. This includes galaxies that have been "ripped apart" by an encounter with 
another galaxy, as well as smaller satellite galaxies like the Magellanic Clouds. [8] 

This classification of galaxies, and also the distinction between Population I and II stars, is 
rather incomplete and coarse. There are more sophisticated classifications that take account 
of more distinct features, e.g. rings in spiral galaxies [23]. The Hubble classification, however, 
is absolutely sufficient to outline the basic features of galaxies. In this chapter, I will be 
concentrating on (barless) spiral galaxies which are observationally interesting because of 
their clearly defined rotating disk. 

3.1.2 Anatomy of spiral galaxies 

Spiral galaxies have a lot more obvious substructure than either elliptical or irregular galaxies. 
While ellipticals appear to be smooth and featureless, spiral galaxies exhibit several visible 
parts: The aforementioned disk which contains the spiral arms, a central bulge, in the cases of 
SB galaxies a central bar, and finally a spheroidal halo that extends beyond the disk. Recent 
observations have shown that generally there is also a supermassive black hole at the center 
of most if not all spiral galaxies. I will use figures of our Milky Way galaxy to illustrate the 
dimensions of a typical spiral galaxy. 

Disk and spiral arms. The differentially rotating disk contains the active star formation 
regions which form the spiral arms. Thus, it contains mostly young stars (Population I), as 
well as gas and dust [8, §1.1]. The spiral arms are not a manifestation of a differentially 
rotating disk, as one might naively think, but are thought to be the result of two coexisting 
mechanisms: quasi-stationary density waves (according to the Lin-Shu hypothesis [54]), and 
self-propagating star formation regions where supernovae explosions create blast waves that 
trigger compression in nearby gas and dust clouds and therefore initiate star formation [61]. 
Usually, the spiral arms trail the rotation of the disk, but in some cases, e.g. NGC 4622, 
the spiral arms are leading, i.e. pointing into the direction of the disk's rotation [61]. 

The measured rotation velocities are constant in the outer regions and even extend be- 
yond the visible edge of the disk, indicating that the mass-to-light ratio T is not constant 
throughout the galaxy [8, §10.1]. A detailed discussion of the observed rotation curves follows 
in §3.2. 
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Figure 3.2: The Sombrero galaxy^ (M104, Type Sa) has a large and extended central bulge 
of stars, and dark prominent dust lanes that appear in a disk that we see nearly edge-on. 
Billions of old stars cause the diffuse glow of the extended central bulge. Close inspection 
of the bulge in the above photograph shows many points of light that are actually globular 
clusters. M104's spectacular dust rings harbor many younger and brighter stars. The very 
center of the Sombrero glows across the electromagnetic spectrum, and is thought to house 
a large black hole. 

The characteristic thickness^ of the disk is different for older and younger stars. While the 
characteristic thickness for younger stars is about 200pc in the Milky Way, for older stars, e.g. 
our sun, it is about 700pc. In general, the thickness of the disk is small compared to its radial 
extent, which is characterized by the Holmberg radius^ of about 17kpc for the Milky Way [8, 
§10.1]. Hence, the disk is usually modelled as a thin (two-dimensional) disk. Observationally 
it was found that the surface brightness / of spiral galaxies follows an exponential law, 

/(r) = /oe-'■/^^ (3.1) 

where the disk scale length = (3.5 it 0.5) kpc for the Milky Way [8, §1.1]. In this plane, 
the disk population follows near-circular orbits around the center of the galaxy [61]. 



Bulge and bar. According to Wyse et al. [94], the bulge is commonly defined by allocating 
all "non-disk" light into the bulge. Generally, it has a nearly spherical shape and the flattening 
is consistent with the slow rotation of the bulge. It consists mainly of old Population II stars 
and is metal-poor compared to the typical disk population. The bulge of the Milky Way has 



^Figure reproduced with permisson from the Hubble Heritage Team (AURA/STScI/NASA). Figure and 
(abridged) text taken from http://antwrp.gsfc.nasa.gov/apod/ap031008.html 

^The characteristic thickness of the galactic disk is defined as the ratio of the disk's surface density to its 
volume density at the galactic plane [8, §1.1]. 

''The Holmberg radius is defined as the radius of the isophote with a surface brightness of 26.5 mag arcsec"^ 
in the B band, which roughly corresponds to l%-2% of the sky's brightness [8, §1.1]. 
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an effective radius^ of rg ~ 2.75kpc [94] and therefore extends above and below the galactic 
disk. 

The term "bar" is widely used in the literature, generally without being defined precisely. 
In most cases the bar is regarded as an elongated bulge from whose end the two main spiral 
arms originate [94]. Others consider the bar to be part of the galactic disk [61]. Independent 
of the precise definition of the bar, it certainly breaks the azimuthal symmetry of galaxies 
near the center, that is otherwise given if one assumes the density perturbations of the spiral 
arms to be reasonably small. 

Central black hole. It has long been suspected that the matter concentration at the center 
of a galaxy is high enough to harbour an extremely massive black hole. The center of the 
Milky Way is marked by the radio source Sgr A*. Indeed, recent observations of several 
stars that orbit Sgr A* have shown that a supermassive black hole is likely to be located at 
the center of the Milky Way. The technically correct deduction from the observations is that 
a mass of (3.7±1.5) X lO'^A^Q is located within a galactocentric sphere of radius 124 AU [78]. 

More recently, the diameter of the radio source Sgr A* has been observed directly and is 
estimated at only ~ 1 AU [84] . This provides very strong evidence that there is an extremely 
compact object at the center of the Milky Way. Such an object is commonly interpreted as 
a supermassive black hole, although there are ongoing discussions regarding whether other 
objects can bind that much mass in such a small volume, see e.g. §4.2 and §4.3. 

Halo. The halo and bulge were long thought to be a single spheroid that follows a surface 
brightness profile oc r~^/^. In 1980, Bahcall and Soneira [3] tried to model the gravitational 
potential of the Milky Way with contributions from the disk and the bulge only. The density 
profiles were assumed to be given by the observed luminosities and a constant mass-to-light 
ratio T for each the disk and bulge. This model did not reproduce the observed "flat" rotation 
curves that exhibit a constant velocity for the largest observable radii. Only after adding an 
unseen spherical halo component, with a density profile that falls off as 

p oc r~2 (3.2) 

for large distances, did their model show the flat rotation curves. Using high-velocity Pop- 
ulation II stars in the halo as an indicator for the gravitational field, the radius and mass 
of the spherical halo were bound to be at least -Rhaio ^ 41kpc and Mhaio ^ 4 x 10^^ [8, 
§2.7]. This shows that the galactic halo extends well beyond the disk and is at least six times 
as heavy as the visible parts of the Milky Way. This is one of the many indications of dark 
matter, a form of matter whose existence is inferred only from its gravitational effects [8, 
§10.1]. 

Within the last 15 years, it became apparent that the brightness of the bulge falls off 
steeper than that of the visible halo, indicating as well that bulge and halo should be consid- 
erent different structures [94]. The halo consists of even older Population II stars than the 
bulge, which form globular clusters and a very thin population of scattered individual stars 
[61]. 

For the last decade and longer, several different methods to measure the mass and extent 
of the Milky Way's dark matter halo have been employed. Most of them are discussed 
in a comprehensive review by Zaritsky (1999) [95]: (i) measurements of the disk rotation 
curve [30], (ii) mass estimates based on the escape speed of the halo [48], (iii) statistical 
analysis of the motion of satellite galaxies and globular clusters [96], (iv) timing arguments 

^The effective radius of the bulge is defined as tlie radius of the isphote containing half of the bulge's total 
luminosity [8, §1.1]. 
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for the dynamical system of the Milky Way and the Andromeda galaxy (M31) [29, 68, 83], 
(v) analysis of orbital parameters of the Large Magellanic Cloud [55] and (vi) motion of 
satellite galaxies of other galaxies similar to the Milky Way [97]. A comparison of these 
measurements is displayed in fig. 3.3 and proves to be consistent with a halo model that 
exhibits a mass distribution which rises linearly with the radius. Although the figure shows 
data for up to ~ 300kpc, the Milky Way's halo cannot extend much beyond 200kpc due to 
the presence of our neighbour, the Andromeda galaxy. This should lead us to the conclusion 
that spiral galaxy halos are roughly spherical on a local scale and generally overlap on scales 
that include neighbouring galaxies [95]. 




Radius (kpc) 



Figure 3.3: A comparison of various measurements of the enclosed mass of the Milky Way 
halo^. The solid lines denote the enclosed mass for dark matter halos which cause a constant 
circular velocities of 180 and 220km/s. The data points are generally given with their 90% 
confidence intervals. The methods in chronological order (see text): Einasto & Lynden-Bell 
(iv, 1982), Zaritsky et al. (iii, 1989, ZOSPA, 3 data points), Fich & Tremaine (i, 1991), 
Zaritsky & White (vi, 1994, 2 data points), Lin et al. (v, 1995), Peebles and Shaya et al. 
(both iv, 1995), Kochanek (ii, 1996, 2 data points). 

While the shape of the halo is generally assumed to be spherical, the actual observational 
situation is currently not conclusive [39] and suggests that spherical as well as slightly oblate 
or prolate halos are possible. Recently, more data on the debris of the Sagittarius dwarf 
satellite galaxy suggests that the dark matter halo might actually be slightly prolate [40] 
while earlier microlensing observations suggest that the halo is more likely to be oblate [75]. 
Despite numerous articles that suggest techniques to measure the shape of the dark matter 
halo and others that predict halo shapes from numerical simulations, actual observations 



Figure reproduced with permisson from Zaritsky (1999) [95]. 
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are very sparse. Hence, the observational situation can be described as somewhat vague, 
although most observations and predictions are consistent with a (at least close to) spherical 
halo. 

Most of the observational evidence that was presented in this section applies only to the 
Milky Way, but the overall picture that has been drawn of a spiral galaxy is consistent with 
observations of other spiral galaxies apart from the Milky Way: The disk with its spiral arms 
and the bulge are clearly visible in other galaxies. Observations of nearby galaxies indicate 
that a central black hole is likely to be a standard feature of spiral galaxies [51]. A large 
sample of rotation curves (~ 1100) showed that the flat rotation curve behaviour is a universal 
feature of spiral galaxies [72] which strongly suggests that a dark matter halo exists in every 
spiral galaxy. The same conclusion can be reached from large samples of weak gravitational 
lensing by galaxies [16]. 

After outlining the structure of spiral galaxies, I shall now go on to discuss the contribu- 
tions of the individual components to the gravitational field. 

3.1.3 Contributions to the gravitational field 

All of the aforementioned components of a galaxy have mass and thus contribute to the 
gravitational field of a galaxy. Therefore, constructing a comprehensive model that respects 
the effects of all parts is rather difficult. The dispute over which density profile to use for the 
individual components adds to the difficulty, and obscures an objective comparison of different 
galaxy models. Fortunately, the contribution of all parts is limited to certain distance scales, 
so that not all components have to be included in a reasonably realistic model. A model of 
the very center of the Milky Way, for example, should include the central black hole and the 
dense central region of the bulge [78] , while the effects of spiral arms or the dark matter halo 
almost certainly play no vital role. 

The gravitational influence of the central black hole is noticeable only within a few parsec 
of the galactic center [78], and can be completely neglected for studying the dynamics in 
the disk or the outer regions of the bulge. The bulge and bar dominate the innermost few 
kiloparsec of the galaxy and depending on their relative size, have to be accounted for in 
models of the galactic disk's dynamics. The spiral arms are density perturbations of the disk 
and thus only influence the gravitational field locally. When studying the motions of stars 
or gas in the disk, one should at least discuss the perturbations of the spiral arms, if they 
are not explicitly included in the model. Gas or dust clumps in the disk can have a similar 
infiuence on local disk dynamics [31]. 

The main contribution to the gravitational field in most of the disk originates, however, 
in the dark matter halo. While the disk mass definitely has to be accounted for, the major 
contribution in the outer region of the galactic disk is due to the halo. Dwarf galaxies are 
almost completely dominated by the dark matter halo and even the disk dynamics of the 
largest spirals cannot be modelled accurately without the halo [72]. 

When modelling the gravitational field beyond the visible disk, out to several hundreds 
of kiloparsec, the dark matter halo is practically the only source of the gravitational field, see 
fig. 3.3. 

Summarizing, I note that the dark matter halo is generally the main contributor to the 
gravitational field of spiral galaxies. For very large galaxies, the disk mass plays an important, 
non-negligible role and the central region of a galaxy cannot be modelled acurately by the 
halo alone. 
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3.1.4 Dark matter halo models in general relativity 

Choosing a model for the gravitational field of a galaxy depends of course on the planned 
application, i.e. the aspect of the research to be undertaken. In this thesis, I want to 
concentrate on the general relativistic aspects of the dark matter halo. 

Strong gravitational fields, that require general relativity to describe the resulting dynam- 
ics accurately, occur only close to massive, compact objects, like stars and black holes. On 
a galactic scale, however, one cannot keep track of all ~ 10^^ individual stars and generally 
introduces a locally averaged density function. Because the absolute value of that density 
is very low^ in galaxies, the induced gravitational fields are also weak. So why use general 
relativity to describe the weak gravitational field in galaxies? 

Newtonian physics is obtained in the limit of general relativity where: 

(i) the gravitational field is weak, 

(ii) the particle speeds involved are slow compared to the speed of light, and 

(iii) the pressures and fluxes and stresses are small compared to the mass-energy density 
[60, §17.4]. 

While there is no doubt that all objects that make up a galaxy satisfy the conditions 
(i) and (ii), at present no one knows much about the nature of dark matter and hence, the 
possibility of dark matter being a high pressure fluid, or some sort of unknown field with high 
field tensions, cannot be excluded. It is this uncertainty about the dark matter's equation of 
state (see §2.6) that makes it necessary to include general relativistic effects. Indeed, I will 
show in the following section that many different dark matter candidates exist which exhibit 
pressures that are comparable to or even dominate the energy density. 

When performing measurements of the gravitational field of a galaxy, condition (ii) might 
play an important role. For rotation curves, the tracer particles surely exhibit subluminal 
speeds, but for gravitational lensing measurements, the particles directly influenced by the 
gravitational field are the photons of a background object which per definition travel at the 
speed of light. I will show how this influences the interpretation of gravitational lensing 
measurements in §3.3. 

If one does not require condition (iii), general relativity does not quite reduce to Newton's 
gravity. It is a standard result that using only condition (i) the gravitational potential $ is 
generated by the tt-component of the Ricci tensor [60, §17.4]: 

= Ru . (3.3) 

Invoking the Einstein field equations (2.33) with a stress-energy tensor of the form (2.62) 
gives Ru- 

Rtt= 4tt{p + pi+P2+P3) , (3.4) 

where the pi are the principal pressures obtained by diagonalizing the space-space part of the 
stress energy tensor. Thus, in the case of a perfect fluid {p = pi = P2 = P3 = Pr = Pt) we find 
the field equation 

V^^ = A7r{p + pr + 2pt) =4Tr{p + 3p), (3.5) 

which only reduces to Newtonian gravity if p <C p, i.e. condition (iii). It is now quite obvious 
from (3.5) that the gravitational field is highly sensitive to the pressure if density and pressure 
are of the same order of magnitude. 

^For example, the density in the solar neighbourhood is po = (0.18 ± 0.03) A4q pc^^ (also called the Oort 
limit). [8, §4.2.1] 
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A different reason to examine galactic dark matter in a general relativistic framework is 
the boundary problem. Since there is no observable boundary of the dark matter halo (see 
fig. 3.3), it is possible that the halo is actually a phenomenon of cosmology. In that case, 
since the halo should then smoothly merge with a FLRW cosmology, a proper description of 
the halo is only possible in general relativity. 

As a first approximation to the dark matter halo, I choose a spherically symmetric space- 
time of the form (2.32). As I already discussed in §3.1.2, most of the available measurements 
and simulations of halos are consistent with a spherical distribution of dark matter. Even if 
the halo was slightly aspherical, the high symmetry of (2.32) simplifies the algebra tremen- 
dously and illustrates the results without obscuring them in too complicated a formalism. 
The basic principles that are discussed in this thesis, however, are basic in their nature and 
not specific to spherical symmetry. 

The currently available data is still too vague to definitely identify a possible ellipticity 
of the halo. If it should turn out eventually that spherical symmetry is not suitable for an 
approximate description of the gravitational field of a galaxy, the results presented in this 
thesis can easily be adapted to a system with less symmetry. 

Before we go on to see how the dark matter halos connect relativistically to the actual 
observations in §3.2 and §3.3, I will present some speculative candidates for dark mass-energy. 



3.1.5 The nature of dark matter: speculative models 

As Binney and Tremaine [8] point out, there is no a priori reason to assume that mass and 
luminosity should be well correlated. Nonetheless it was long believed that in the universe, 
there is only matter where there is light. The first indication that there is more matter than 
just the visible kind, was found by Zwicky in 1933 [98]. He was observing the dynamics 
of galaxies in the Coma cluster and came to the conclusion that there is some 400 times 
more mass in that galaxy cluster than indicated by the luminosity of the galaxies and their 
corresponding mass-to-light ratio alone. 

Since then, many more different indicators for dark matter have been found, see [8] for 
a comprehensive overview. The question "What is dark matter?" , however, remained unan- 
swered. Over five decades after Zwicky's discovery, Binney and Tremaine [8] considered the 
nature of dark matter the "probably single most important unresolved question in extra- 
galactic astronomy". Today, we still don't know the answer. 

Many different models exist for dark matter, all of which share the only known property: 
dark matter generates gravitational fields and does not emit photons in any measurable 
quantity. I am going to present the currently most discussed dark matter candidates and 
focus on their equation of state and typical pressure to density ratio, 

Pr + 2pt p Y^iPi rna\ 

w = — = - = — — . (3.6) 

3p p 3/9 

Please note that in this definition, w is not necessarily a constant, but depending on the 
form of matter at hand, can vary with density and pressure. From equation (3.5) we see that 
general relativistic effects can only be neglected in observations if \w\ <C 1. In §3.2 and §3.3 
I will argue that a significant w can in principle be observed when one keeps track of the 
pressure terms in a general relativistic analysis of the data. 



Current estimates of cosmological density parameters. Although the composition of 
dark matter is not known, reasonably tight constraints on the amount of different matter types 
exist today. A major boost in accuracy was achieved by fitting the cosmological standard 
model to the data of the Wilkinson Microwave Anisotropy Probe (WMAP) in 2003 [6]. The 
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density parameter for a matter species x is given in terms of the critical density , = 
(px) //Ocrit, where {px) indicates the average over the whole universe. The first Friedmann 
equation can then be written as 

ntot = Y,^^ = 1 + ^772, (3.7) 

where the sum runs over all species, ao is the curvature radius^ of the Priedmann-Lemaitre- 
Robertson- Walker metric, Hq is the Hubble parameter^ and k G {—1, 0, 1} determines whether 
the FLRW space is a hyperbolic 3-plane, Euclidean or a hypersphere. The Hubble parameter 
is often given in terms of h, which is defined by Hq = 100/ikms~^ Mpc~^. 



Parameter 



Symbol Value 



Total average density of the universe fitot 

Density of cosmological constant A Q\ 
Total matter density 

Baryon density Ob 

Non-baryonic matter density ^^nbm 
Cosmic density of luminous matter (§20.3) ^^lum 

Neutrino density (§21.1.2) Jljy 

Radiation density (§21.1.4) ^^rad 



1.02 ±0.02 
0.77 ±0.06 
0.25 ± 0.04 
0.043 ± 0.006 
0.21 ±0.04 
0.0033 ± 0.0002 
~ 0.001 

0.000046 ± 0.000004 



Table 3.1: The cosmological density parameters as given by the Particle Data Group (2004) 
[67]. Fitting the first acoustic peak of the cosmic microwave background (CMB) angular 
power spectrum yields the values of Jltot) ^^m and fib (§21.1.4 and §21.3.3 of [67]). The dark 
energy density JIa is derived from f^tot — ^m, and 0,^hm is given by — (§22.1.1). All 
values and cr-confidence intervals were obtained using a value of /i = 0.73 ± 0.03. 

The Particle Data Group [67] bi-annually provides the current best values of the density 
parameters, which are listed in table 3.1. The data suggest that the most mass-energy in the 
universe is present in the form of dark energy, which is most typically modelled as cosmological 
constant with w = —1 [67, §21.2.3]. The majority of matter in the universe is dark since the 
luminous fraction of matter is only r^ium/^m ~ 0.013. Although the baryonic dark matter is 
dominating the visible matter, r^b/riium ~ 13, most of the dark matter must be present in 
non-baryonic form: r^nbrn/^m ~ 0.83. 

So, not only do we not know what dark matter actually is, it also looks like there are at 
least two different kinds of dark matter - if the cosmological standard model is correct. And 
as of 2004, it fits the data remarkably well [67, §21.4]. I will now present some suggestions 
what dark matter actually might be. 



Baryonic dark matter could be present as unobserved low-luminosity stars that do not 
have enough mass to burn hydrogen (brown dwarfs, M < 0.08 A^o) or stellar remnants of 
high mass stars, that now would be present as white dwarfs, neutron stars or black holes. 
Since all these objects do emit a small amount of in principle detectable radiation, their 
density can not exceed l^steLrem. ~ 0.03 [8, §10. 4. la]. 

*The critical density, which is also known under the name of the Hubble density, is defined as pcrit = ^J^q ■ 
If the average density of the universe equals pcrit, i.e. f2tot = 1, (3.7) predicts that the universe is spatially 
flat. [8, §9.A] 

^ao and Hq refer to the current epoch. 
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Other examples of massive objects that are very difficult to detect are primordial black 
holes that were created during the big bang, or Jupiter mass objects (~ 10-3 A^o) [61]. Ah 
of the aforementioned objects are generally referred to as MAssive Compact Halo Objects 
(MACHOs) , and detection of their presence and estimation of their total mass is attempted by 
microlensing surveys [66]. Together with MACHOs, cold gas clouds have also been suggested 
as baryonic dark matter candidates [22]. 

All of these baryonic dark matter candidates have non-relativistic speeds and are therefore 
part of the class of Cold Dark Matter (CDM). The parameter w for CDM can be estimated 
by considering an ideal gas of particles with mass nip: 



where is Boltzmann's constant, T is the "temperature" of an ideal gas as given by 



and c the speed of light. Therefore, the non-relativistic particle speed of CDM indicates 
negligible pressure, i.e. w ^ 0. 

Exotic particles. Massive neutrinos as well as Hot and Warm Dark Matter (HDM, WDM) 
are generally ruled out as major dark matter candidates, since their speed is usually too high 
to allow for the dark matter clumping which is observed in galaxies [8, 63]. Consequently, 
if they are to be significant components of the dark matter, exotic particles should classify 
as CDM. The two most prominent exotic CDM particles in the literature at the moment are 
axions and Weakly Interacting Massive Particles (WIMPs) [67, §22.1.2]. 

Axions were first postulated as a means to solve the charge-parity problem of Quantum 
ChromoDynamics (QCD), but they also occur naturally in superstring models [67, §22.1.2]. 
Two experiments that search for axions are currently collecting data, one at the Lawrence 
Livermore National Laboratory, California and the CARRACK experiment in Kyoto. Both 
of them have so far only excluded axions with certain rest masses, and have not yet actually 
found positive evidence for axions [67, §22.2.2]. The current astrophysical observations, and 
laboratory experiments, indicate an upper limit on the axion mass of rua ^ 10~^ eV [67, 



WIMPs are characterized as particles with a rest mass of roughly 10 GeV to a few TeV 
and a weak annihilation cross section that makes them difficult to detect. The currently 
favoured WIMP candidate is a Lightest SuperParticle (LSP) in super symmetric models, and 
amongst the LSPs, the "neutralino" seems to be the most promising one [67, §22.1.2]. With 
both, axions and WIMPs, being candidates for Cold Dark Matter, the expected pressure 
content of dark matter in these cases is negligible compared to the mass-energy density and 
hence, w ^ 0. 

Non-CDM fields. Apart from the currently favoured Cold Dark Matter models, a few 
approaches consider scalar fields and Bose-Einstein condensates that exhibit a significant 
field pressure or tension. Amongst these proposals for DM candidates is a "string fluid" 
which suggests w ~ —1/3 [87] (though other models of string fluids also appear in the 
literature), and many different approaches to Scalar Field Dark Matter (SFDM) which is 
usually modelled through at least one scalar held in a general relativistic framework with 
a self-interaction potential. All SFDM models exhibit signiflcant pressures, e.g. Schunck's 
model has a position dependent w S [0, 1] [80], Peebles considers a scalar field with a purely 




(3.8) 




(3.9) 



p.393]. 
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quartic self-interaction potential and finds w = 1/3 [69], tlie scalar field model of Arbey et 
al. (which is variant of a Bose-Einstein condensate) has a quadratic equation of state and 
thus^*^ w = Xh^ /{A'm'^c)p [1] and Matos et al. claim that their SFDM model even exhibits 
w ^ 10^ [57]! 

Modified physics. A completely different approach to explain the dark matter problem is 
to demand that the theory of gravity must be modified for distance scales as big as those of 
the dark matter phenomena. This idea, generally referred to as MOdified Newton Dynamics 
(MOND), is formulated by Milgrom [59]. Until this theory is extended into a "clean" fully 
relativistic invariant theory, it can not be considered a complete solution to the dark matter 
problem. 

Since there are competing dark matter candidates with different values for w, any way 
of actually observing w would mean that it is possible to further restrict the possible dark 
matter candidates. In the following two sections, I will outline a way of in principle measuring 
w in spiral galaxies. 

3.2 Rotation curves 

The term "rotation curve" generally refers to a set of circular rotation velocities of tracer 
particles in the disk of spiral galaxies that were obtained at different radii. Given the circular 
rotation velocity, one can infer the mass inside a sphere with the corresponding radius. A 
rotation curve then provides information about the radial matter distribution. If one addi- 
tionally assumes a model for the 3D shape of the matter distribution in the spiral galaxy, i.e. 
respecting the different non-spherical contributors, the rotation curve determines the entire 
matter distribution. 

The Doppler shift of emission lines makes the observation of the tracer particles' mo- 
tion possible. Mostly the 21 cm emission line of neutral hydrogen (Hi) is observed at radio 
wavelengths. Hi gas observations generally extend to larger radii than other tracers. [8, 
§10.1.6] 

3.2.1 General shape of rotation curve profiles 

In the early days of rotation curve measurements in the 1950s, the rotation curves were 
modeled after the exponential disk model of galaxies, see (3.1). This model implied that a 
rotation curve would exhibit three characteristic regions: (i) the central region where the 
velocity rises linearly with distance, (ii) a region where the velocity reaches its maximum, the 
so called turnover radius, and (iii) the outer region, which was expected to show a Keplerian 
falloff oc r~^/^ since the galaxy supposedly resembled a point mass at large distances. [8] 

With this expectation in mind, it is not a big surprise that the first observers assumed 
the Keplerian falloff to be there, even though the observation extended only up to a flat 
region of constant velocity, which was interpreted as the turnover radius. It was only in 
the 1970s that the improved instrumental sensitivity allowed observations out to larger radii 
and it was found that the flat region of the rotation curve extended further than agreement 
with the exponential disk model would allow. Instead, the flat region extended as far as 

^"A is the quartic coupling parameter of the scalar field, h is Planck's constant divided by 2tv, m is the mass 
associated with the scalar field, i.e. the boson mass of the Bose-Einstein condensate, and c is the speed of light. 
In SI units, to ~ ( TcF^) (''n?")^ ( 9 sxio-^'^ kg/m^ ) ■ ^^'^ comparison, the density of the solar neighbourhood 
is pQ = (0.18 ±0.03) Ale pc"^ ^ 1.2 x 10"^°^. 
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the observations of the rotation curve reached, without any indication of dechning rotational 
velocity. Consequently, Freeman suggested in 1970 that "there must be [...] additional matter 
which is undetected, either optically or at 21cm [...], and its distribution must be quite 
different from the exponential distribution" [32]. 

Today, several thousand rotation curves have been observed and are publicly available 
in electronic databases [86]. Data analysis of about 1100 rotation curves shows common 
properties for all sizes of spiral galaxies: The most massive galaxies show a slightly declining, 
though non-Keplerian, profile in the outermost region; intermediate galaxies have a nearly 
flat region in the outer parts; and dwarf galaxies show a steadily increasing velocity profile 
throughout the observable region [72]. This analysis brought forth a completely empirical 
universal rotation curve for spiral galaxies that is only parameterized by the total luminosity 
L: 



/ 1 07^1.22 ^2 

with 

a(A) ^ 200^fix AO-4i{o.80 + 0.491ogioA + ^^ffS^}~'^' (3.11) 

6(A) = 0.72 + 0.441ogio A (3.12) 
c(A) = 1.6 e-°-^^ (3.13) 

X = r/Ropt^^^ (3.14) 

V A 13kpc 

A = L/L, (3.15) 

= 1.0 X 1O^°/i-2L0 , (3.16) 

where Lq is the Sun's luminosity in the visual band [8], and h is the dimensionless Hubble 
parameter. The 6(A)-term of the velocity arises from the stellar disk which is modeled as 
an approximation to the exponential disk model (3.1). The dark halo contribution of this 
universal rotation curve is given by the c(A)-term and has an asymptotically constant rotation 
velocity of [72] 

Voo^l.3a{X)e-°-^^ . (3.17) 

This corresponds to a halo density that falls off as p oc and consequently implies an infinite 
total halo mass. Such a model is usually referred to as a "singular isothermal sphere" (SIS). 
The name arises from the solution of the Newtonian equation of hydrostatic equilibrium. 



for an isothermal ideal gas with 



dp ^ M{r)p{r) 
dr ' 

dp ks T dp 



(3.19) 

Combining both equations, one finds using M'(r) = 47r r2p(r): 



dr nip dr 



d / 2 d In p A 47r nip 2 



r 



dr \ dr J ksT 

with the solution 



r>(r), (3.20) 



P{r) = ^-,. (3.21) 
This density profile diverges for r = 0, hence the name singular isothermal sphere. Also the 
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total mass of the SIS diverges linearly with radius. 

Alternatively, numerical simulation of the evolution of CDM halos in a cosmological con- 
text suggest a different density profile, one that falls off as p oc for very large distances. 
This still corresponds to a logarithmically divergent total mass. Such halo models are usually 
referred to as NFW-halos, named after Navarro, Frenk and White who conducted the rele- 
vant simulations [62, 63]. The associated rotation curve is asymptotically declining, although 
the decline only occurs in a region that lies beyond the region in which rotation curves of 
individual galaxies can be observed (> 100 kpc). However, statistical analysis of large sam- 
ples of satellite galaxies in orbit around their primary galaxies seem to support the notion of 
NFW-halos and the declining rotation curve at large distances [73]. Statistical methods that 
use weak gravitational lensing data (§3.3.3) also seem to favour NFW halos [34, 45]. 

Summarizing, the direct evidence of individually observed rotation curves exhibits a gen- 
eral behaviour of an asymptotically flat rotation curve up to the last observable data point. 
This corresponds to a density that falls off as p oc r~^. However, there are indirect indications 
from statistical methods that the rotation curve ultimately declines, which implies a density 
profile that goes asymptotically as p oc r^"^ or even steeper. If the total mass of the galactic 
halo is to be finite, then the density must asymptotically fall off faster than p oc r~^. 



3.2.2 Measuring the rotation curve through direct observations 

To calculate the relationship between the gravitational potential $ and the observed redshift 
Z of the tracer particles, I make two assumptions that will simplify the relevant algebra 
tremendously: 

1. The tracer particles are in circular orbits around the center of the galaxy and are 
confined to the plane of the disk. 

2. The gravitational field is spherically symmetric. 

The first assumption agrees with what is generally known about spiral galaxies (see §3.1.2) and 
thus, the rigid restriction to motion within the plane (6 = tt/2) and circular orbits (dr = 0) 
is acceptable as a first approximation. Assumption 2 is motivated by the approximately 
spherical shape of the dark halo. There are, however, also aspherical contributors to the 
gravitational field, which is in the region of interest mainly due to the stellar disk. But since 
we are for now only interested in the motion of the tracer particles within the galactic plane, 
the symmetry in the plane is circular, no matter whether the gravitational field is spherically 
or only azimuthally symmetric. Hence, I will utilize a metric of the form (2.28) or equivalently 
(2.32). 

The following derivation of the relation between the observable redshift z± and the metric 
function $(r) is based on a modification of the presentation of Lake [52] and Weinberg [92]. 



Constants of motion 

The geodesic equations (2.8) of the metric (2.28) with an affine parameter x are 

d^t A'(r) dt dr , , 

dx A{r) dx dx 

^ dV A'jr) fdty B'jr) / drV r ( dOV rsm^9 (dipV 

dx'^ 2Bir){dxJ ^2B{r)\dxJ B{r) \dxJ B{r) \dx) ' 



d'^e 2dedr . , nfd^V 

-r^ H — sm fci COS 6^ — 

dx '^dxdx Vdxy 



= — + - — —- sine cos9[-^] (3.24) 
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dV 2 dip dr ^ dip d9 , , 

= _| + --^_ + 2cot0-f — , 3.25 
dx r dx dx dx dx 

where the shorthand symbol ' = d/dr has been used. Due to assumption 1, (3.24) is imme- 
diately satisfied, and we can disregard in the following discussion. The remaining geodesic 
equations simplify to 

dh A'(r) dt dr ,^ , 

= TT^ + ^rrr:^:^ 3.26 

dx^ A{r) dx dx 

^ dV / d^y /dr 

dx^ ^ 2B{r) \dx) 2B{r) \dxj Bir) \dxj ^ " ^ 

d'^ip 2 dw dr , , 
dx r dx dx 

Dividing (3.26) and (3.28) by dt/dx and dp/dx respectively, yields 

±fln(A{r)^]] = (3.29) 



dx I V dx 

:T^^ln(r^^')^ = 0, (3.30) 



dx I V dxy 

which motivates the definition of two constants of motion: 

, . ,3.31) 

' . ...(-)|^^^|. (3.32) 

where I stated 7 also in terms of the metric function $(r) of the geometry (2.32). These 
quantities are proportional to the conjugate momenta of the cyclical coordinates t and p, 
where the term "cyclical" refers to their absence in the metric components ga/s- 

With these definitions and multiplication by 2i?(r)dr/dx, the remaining geodesic equa- 
tion (3.27) simplifies to 

which gives yet another constant of motion. 

By inserting all four constants of motion, 9, 7, i and /C, into the metric (2.28), I find the 
proper time r to be given by 

dr2 = -ds2 = /C2 dx^ (3.35) 

This was expected, since proper time is in general an affine parameter for timelike particles. 
The constant /C is not a dynamical constant of motion, but a mere manifestation of gauge- 
freedom to choose the affine parameter. 



The special case /C = =^ dr = is characteristic of the motion of null-particles, i.e. 
photons or other particles that travel at the speed of light. For these particles, (3.34) simplifies 
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to 

= -B{r)( — ] +-L^-_, (3.36) 



dxj A{r) 



r 



where the affine parameter (2.25), Xi is generally chosen for convenience to satisfy 

> 1 for r — > oo , (3.37) 

if the particle actually reaches r = oo at any point of its trajectory. 

All other values for fC invoke a non-vanishing proper time interval dr 7^ and thus 
correspond to the motion of timelike particles. Defining 

7 = /C7; i = lC£, (3.38) 

for all timelike particles, makes (3.34) independent of K, (as long as /C 7^ 0): 

At the same time, the constants of motion 7 and (. gain a familiar interpretation: 

7 = e2*(^) ^ (3.40) 
dr 

is easily identified as the energy per unit mass at infinity and 

i = r'^ (3.41) 
dr 

as the angular momentum per unit mass [60, §25.3]. From now on, I choose fC = 1 for timelike 
particles so that dr = d^. For null trajectories, however, the affine parameter x must still 
be specified. 



Circular orbits of timelike particles 

For circular orbits, r = re is another constant of motion, where r^ labels the radius of the 
emitting particle's orbit. Thus, dr = and using (3.40) and (3.41), the geodesic equation 
(3.27) simplifies to 

where again, I used the metric function <I>(r) instead of A{r). Combining the geodesic equation 
(3.42) with (3.39) for circular orbits gives the energy and angular momentum for the emitter 
at radius r^: 

Y^l - re $'(re) 

4 = ^^P= (3.44) 

Y^l - re $'(re) 

Hence > 0, from the condition 

0<r^>'(r)<l (3.45) 
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for the existence of circular orbits. 



Observing light emitted by particles on circular orbits 

The redshift of light is given by [79, p. 47ff] 



l + z 



(3.46) 



where A is the wavelength, e denotes the event of emission, o the event of observation, 
V"' = dx°^ /At is tangent to the worldline of the emitting particle / observer and = dx^/dx 
tangent to the null geodesic of the emitted / observed light. Furthermore, r is proper time, x 
is a suitable affine parameter for null curves, and Qajj is the contravariant metric. The affine 
parameter for the photon trajectory in a static spacetime (2.25) yields a constant energy at 
infinity (3.31) for the null curve: 



-9tt 



dt = e2*('-) dt 



In 



1 



(3.47) 



Let us assume the observer to be at rest at infinity in an asymptotically flat space time with 
$(oo) = 0, so that the only non-zero component of the observer's tangent vector is 



dx 



dr 



-^(oo) 



(3.48) 



The emitting particle performs a circular orbit, thus dr = d^ = and from (3.41) and the 
invariant proper time interval, dr^ = — ds^ (2.32), follow the non-zero components of the 
emitter's tangent vector 



dx* 
d7 



The observed redshift then takes the form 

Qtt 



l + z 



ri 



ie dip 



(3.49) 
(3.50) 



n 1 dt 



(3.51) 



Inserting the affine parameter dx = —gtt dx* yields 



l + z 



1 + 



9tt rl dt 



N 



(3.52) 



where d(/7/dt|^ refers to the null path of the light, conveying information about the emitting 
particle to the observer. Recalling the constants of motion 7jv and In given by (3.31) and 
(3.32), which are defined along the null trajectory of the photon, I define a new constant of 
motion, the impact parameter 
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The indices N emphasize the fact that b is constant for all points along the null worldline of 
any photon. 




Figure 3.4: A particle with mass m on a straight trajectory with velocity v and position 
vector r has the impact parameter b. The absolute value of the particle's angular momentum 
is given hy L = Imf x v\ = mvr sin (p. 



To see that b is the impact parameter at infinity (i.e. the observer), remember the classical 
definition of the impact parameter (fig. 3.4) 



L = Imf X v\ = mvr sin(<(r, v)) = mvb 



b 



L 



mv p 7 



(3.54) 



Please note that for photons j = pc and in geometrical units c = 1. Since it is not possible 
to measure r^ directly, it is necessary to rely on observing the apparent radius, which is the 
impact parameter b at infinity. 

It follows now from (3.52), using (3.44) for the angular momentum of the emitting particle, 
ie, and the metric components g^^ip and gu of (2.32) at the radius r^ that 



1 + z 



1 + 



r|<l>'(re 



e-1>(r-e) + 



r2(l -re$'(re)) 



r-g rey/re^'{re) dcp 
-emre) rl^- re$'(re) 



dv9 



N 



^1 - re$'(re) re Vl - r.^'ir,) At 



N 



In terms of the impact parameter b, given by (3.53) for the event of emission (r^r = r^), this 
is finally 

l + z±= \,, . 1 - ^'"'^•'^^ 1 • (3.55) 



1 ±\b\y^re^'{re) 



a/1 - re$'(re) \ e 



The terms involved in (3.55) have a simple interpretation: z± is the total redshift of light 
that was emitted by particles moving on circular orbits in a geometry of the form (2.32), 
where z^ is the redshift of approaching and z- that of receding particles. The second term 
vanishes when the light ray moves radially outwards {dx'^/dx = 0), perpendicular to the 
circular trajectory of the emitter, which can easily been seen from (3.51). Therefore, it 
represents Doppler redshift. The first term is independent of the emitter's motion and only 
dependent on the the gravitational potential $(re). It represents gravitational redshift. Thus, 
I define gravitational redshift 
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and Doppler redshift (where b > for the rest of the discussion) 

2 re^l-re^>'(re) ' 



(3.57) 



Similar derivations of Zg and Z can be found in [52, 64]. 

Through (3.40), (3.41), (3.43) and (3.44), these terms are also related to dt/drlg and 



1 + Za 



Z 



7ee 



~2<I>(r,) 



dt 
d7 



dip 
d7 



(3.58) 
(3.59) 



From the proper time interval for the worldline of the emitter's circular orbit, we get an 
expression for dt/drl^: 



dr^ = e2*(^^)dt2-r2dv?2 
dt ' 



(3.60) 



dr 



1 + ri 



^4 



dip 

dip 
d7 



+ 



difV 



(3.61) 



Since the first term of the particle's "circular velocity", rg dip/drl^ <C c, occurs only in 
the second order in (3.58), one can already suspect that the gravitational redshift will be 
negligible compared to the Doppler redshift (3.59) in which the "circular velocity" occurs in 
first order. The quotation marks are necessary, as at the moment the different quantities b 
and re are compared. To get a valid comparison, their relationship has to be examined. 



Impact parameter b vs. coordinate radius r 

For legibility the index e is now dropped and all radii r refer to the radius of the emitter's 
circular orbit. 

Generally it is not possible to measure Z(r) directly. This is firstly due to the fact that the 
observer measures the position of a signal in terms of the apparent radius or impact parameter 
b. Secondly only the total redshift z±{b) is observable. Finally, the observer measures many 
different redshift values from various emitters along the line of sight, given by b. A mapping 
r{b) for a corresponding redshift value z±{b) is necessary to filter out the appropriate z±{b) 
from all available values. Thus, a way to obtain the gravitational potential $(r) is 

1. Measure total redshifts z^(b) and z^(b). 

2. From that calculate Z{b) = ^ {z^{b) — z+(6)). 

3. Obtain <^{b) from Z{b). 

4. Interchange variables b <^ r to get <l>(r). 

This procedure is of course idealized and in practice, one faces problems that arise from noisy 
data and imperfect symmetry between the approaching and receding emitters for the same 
impact parameter, z±{b), which introduces even more noise in step 2. Other problems that 
arise due to insufficient modeling are discussed in §3.2.4. 
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In the given scenario, the observer is looking at redshifts of an edge-on galaxy. It is 
therefore not obvious which radius belongs to a redshift value that originated from some 
emitter along the line of sight. Hence, a way of relating a characteristic redshift to the 
emitter's radius is needed. 

I will now establish the relationship between the coordinate radius and the impact pa- 
rameter, r(b), which is needed to perform step 4. This will also identify which redshift value 
z±{b) should be picked. To find the corresponding radius for the redshift values of a given 
impact parameter b, I square (3.57) to get 



?1 

62 



r(l - r^' 



(3.62) 



Now assume this quantity to be monotonically decreasing with increasing r, which is synony- 
mous with 

> r$" + 2r($')^ • (3-63) 

The maximum value of Z'^ for a fixed impact parameter b is then related to the minimum 
value of the radius of the observed null geodesic, r = = rAr^min ^ dr^v = 0, see fig. 3.5. 




dr = 



Figure 3.5: This sketch shows the trajectory (solid line) of all photons that are observed with 
an impact parameter b. The dash-dotted line marks the line of sight corresponding to b. The 
observed photons could have been emitted from any point on the solid line. However, from 
(3.62), one can conclude that photons emitted at the point of closest approach, rjv min) to the 
coordinate origin exhibit the strongest Doppler shift. See also fig. 3.6. 



It then follows from the lightlike line element 



that 







9tt - 


d(/? 




>/"'" 




N,rmi\ 


V 9w 
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A'',min 



P^Cf-Ar.min) 



(3.64) 



(3.65) 



Using the definition of the impact parameter (3.53), one obtains the mapping r[b) from (3.65) 
for the moment of emission at the radius r = ri\j ^i^: 



dr 



(3.66) 
(3.67) 



Since (3.62) is assumed to be monotonically decreasing, (3.66) and (3.67) are only valid for 
the maximum redshift value of Z for a given impact parameter b. 

To perform step 3, I want to find a relationship only involving Z{b), b and $(6). Hence, 
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eliminate r in (3.62), using (3.66) and (3.67), so that after careful substitution the result 

Z^ = .i£e-»)=4(-S^l (3.68) 



db 

is obtained, which corresponds nicely to the analogous relation in the Newtonian limit^^, 

Z =rN- (3.69) 

drAT 

where rjy ~ & is the observed radius and e"^*^^^ ~ 1. 

In principle, the potential ^{b) can now be obtained from the measured quantities Z and b, 
by integrating using an appropriate boundary condition for ^{b). Once ^{b) is determined, 
$(r) is given by the mapping (3.66) between b <^ r. From here on, I assume that $(r), 
and especially $'(r) = d$/dr, has been extracted from the data and can be considered as 
"known" . 

Since the relation between r and b is now available, it is possible to compare the gravita- 
tional and the Doppler redshift directly. First I define the circular velocity Vc = [rd(/7/dr]g 
and then apply the mapping (3.66) to (3.59): 

1 + z, ^ e-*w|l + i^,2 + OK]} (3.70) 

Z = e-*(")t;e. (3.71) 

For e*^'^) « 1 it is now clear that Zg <^ Z since the rotation velocities have typical values 
around Vc ~ 250km/s <C c = 1 (in geometrical units) and thus, the gravitational redshift can 
be neglected. 



3.2.3 Observational techniques 

Astrophysical observations are usually carried out without considering general relativity and 
as I have shown with (3.71), this approximation is very good to first order in the rotation 
velocity Vc- The derivation outlined in the previous section, however, considers only edge-on 
galaxies, which posit only a fraction of the observed spiral galaxies. Due to random orientation 
in space, most spiral galaxies are somewhat inclined to the line of sight. I will now outline 
actual observational methods, starting with nearly edge-on observations, followed by methods 
of observing the rotation curve for the more common case of inclined spiral galaxies. 

All direct observations of rotation curves use the Doppler shift of a well known emission 
or absorption line to determine the emitting particle's velocity component along the line 
of sight. Typical lines that are utilized are those of Ha, CO, especially neutral hydrogen 
Hi, and other optical emission lines [86]. Since there are many emitting particles along 
any line of sight, and different emitters generally have slightly different velocities, a whole 
redshift- / velocity-profile is obtained for each observed position. The different methods of 
data reduction attempt to estimate an appropriate average rotation velocity. 

For nearly edge-on galaxies, the Envelope- Tracing method, also called the Terminal- 
Velocity method, is used. As illustrated in fig. 3.6, it establishes the maximum or terminal 
velocity along the major axis of the observed galaxy. To do that, the maximum intensity /max 
and the noise level /noise are determined from the velocity profile, and from that a threshold 
intensity /threshold is calculated which corresponds to about 20% of the maximum intensity. 



^'^The index A'^ refers here to the Newtonian limit, not the null curve. 
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'terminal 



A 



Figure 3.6: A wavelength- (A) (or analoguously velocity- or also redshift-) profile showing the 
signal intensity I for a specific impact parameter b. The redshift ^terminal is determined to 
correspond to the radius of closest approach for b. The other symbols are explained in the 
text. 

The velocity that occurs with the threshold intensity near the edge of the velocity-profile 
is defined to be the terminal velocity .^terminal and then corrected for systemic velocity, in- 
clination, velocity dispersion of the observed gas and instrumental resolution to obtain the 
rotation curve [86] . This is essentially the same as picking the maximum velocity as required 
by the derivation in §3.2.2 plus observationally necessary corrections. See also fig. 3.5. 

In galaxies with a reasonable inclination, the line of sight does not run through the 
plane of the disk, but pierces it from top to bottom. The smeared out distribution in the 
velocity-profile for any observed position then originates from the thickness of the disk and 
the random and anomalous velocities of the emitters within, see fig. 3.7. A representative 
rotation velocity is generally obtained by an intensity-weighted velocity average of the entire 
profile. In the outer parts of the galactic disk, the line profile can be assumed to be rather 
symmetric around the peak-intensity value and therefore, the intensity-weighted velocity can 
be approximated by the velocity with the highest intensity [86]. 



Figure 3.7: A galaxy (grey ellipse) which is inclined to the line of sight (solid) is shown from 
the side to illustrate the finite thickness effects of the disk. The wavelength/velocity profile 
(below) arises from the Doppler shifts of the different velocities of the emitters along that 
line of sight. The emitters which are closer to the center of the galaxy (right dotted line) 
have a higher circular velocity. Therefore, in case they are receding from the observer, they 
appear redder. 




A 



To increase the reliability of the obtained rotation velocities, the initial rotational curve, 



38 



CHAPTER 3. GALACTIC DARK MATTER HALOS 



RCO, can be used to reconstruct the measurable velocity profiles. The difference between 
these profiles and the actually observed ones yields a corrected rotation curve RCl. This 
procedure can be repeated until the difference between the observed and corrected velocity- 
profiles becomes minimal and stable [86]. 

If the rotation velocity is obtained for the whole surface of an inclined galaxy, one can 
apply a tilted ring fitting procedure to obtain a rotation curve that relates only radius and 
circular velocity [31]. This averages all obtained rotation velocities of a certain radius and 
maximises the amount of data that was used to obtain the rotation curve. 

Since the actual observations use more data than that which can be obtained using the 
method in §3.2.2, I will now go on to extend the general relativistic model to account for 
inclined galaxies and systemic velocity. 

3.2.4 Observationally necessary improvements of the model 

In §3.2.2, I presented the redshift of light emitted by a source in the geometry (2.32), as it is 
detected by an observer resting in the galactic plane {9 = tt/2) at infinity. 

Real life galaxies usually show a radial systemic velocity away from the observer on Earth 
(or at least the solar system) . This is generally due to the Hubble expansion of the universe 
and peculiar motion of the observed galaxy, both of which lead to an additional shift of the 
observed emission/absorption line. 

Most observed galaxies are also inclined to the line of sight and therefore the observer is 
not within the galactic plane. This affects the Doppler redshift of the observed light. 

In this section, I am going to introduce these further conditions to the total redshift 
(3.55). Finally I will argue that it is generally not necessary to consider general relativistic 
effects when obtaining the Doppler redshift from galaxy observations. 

Inclination correction 

To see how the inclination angle affects the Doppler redshift, I separate the observer and 
galaxy coordinate system. The galaxy coordinate system (x, y, z) is still aligned with the 
galactic disc and the coordinates of any particle therein shall be given by 

X = f cos (p (3.72) 
y = r sin 
z = 0. 

In his or her coordinate system, the observer is thought to sit on the x-axis at infinity, i.e. 
(r = oo, 9 = 7r/2, ip = 0). Thus, the observer coordinate system (x, y, z) has the same origin 
as the galaxy coordinate system, but the relative orientation between the two depends on the 
galaxy's inclination angle ip and the position of the observed particle. 

To find the appropriate transformation rules between both coordinate systems, two rota- 
tions have to be performed: The first rotation about the y=y-axis accounts for the observed 
inclination angle V'j see fig. 3.8. The second rotation is necessary since for convenience, I 
want the connecting light ray between the emission event and the observation to lie in the 
observer's 9 = -7r/2 plane. This will be a rotation about the observer x-axis (fig. 3.9) that 
makes sure that the observed particle's z-coordinate vanishes, i.e. 9 = 7r/2. 

The inclination by an angle tp between the x- and x-axis is represented by a simple rotation 
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Figure 3.8: The left image shows the galaxy coordinate system {x, y, z) as viewed from the 
y-axis. The galaxy (grey ellipse) is inclined to the line of sight (dash-dotted line) with an 
angle ip. The right image shows the galaxy with its coordinate system (dashed) in the new 
intermediate coordinate system (x, y, z) as viewed from the y-axis. The line of sight now 
coincides with the x-axis. 
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about the coinciding y=y-axis (fig. 3.8): 

cos ^ COS 

sin(^ I . (3.73) 
— sin cos (p 

The observer x-axis is the line of sight to the center of the observed galaxy. Hence, a rotation 
about this axis does not change any results obtained by that observation and can be used to 
simplify the geometric constraints. Such a rotation (fig. 3.9) by an arbitrary angle C is given 
by 

cos (p cos V' \ 
sin COS C + COS sin sin -0 . (3.74) 
sin (p sin — cos p) cos C, sin ^\} j 

Both angles, '0 and are Euler angles due the rotation about the y- and x-axis. The third 
Euler angle is zero, as a rotation about the z-axis is not necessary. Since the ^;-component 
has to vanish, the rotation angle C is given by 

= sin (p sin (" — cos p cos ^ sin 
sin -0 

^ tanC = -. 3.75 

tan 99 

So C is just the angle between the apparent major axis of the observed galaxy and the 
connecting line between the galactic center and the observed particle. Finally, a particle in 
the galactic plane with coordinates (3.72) has the following spherical polar coordinates in the 
observer system: 

r = f (3.76) 



fy\ ( tan-f/^X 

p = arctan I — I = arctan 

V X / \ sin ( J 
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Figure 3.9: The left image shows the inchned galaxy in the intermediate coordinate system 
{x, y, z) as seen by the observer. The small solid circle marks the position for which the red- 
shift is to be determined. The right image shows the galaxy with the intermediate coordinate 
system (dashed) in the final observer coordinate system (x, y, z) as seen from the x-axis. The 
galaxy has been rotated by an angle Q so that the observed emitter (circle) lies on the y-axis. 



where C is given by (p and tjj through (3.75). 



Corrected redshift relation 



Now I apply the additional conditions of inclination -0 and systemic radial velocity Vg to the 
redshift relation (3.51). In the observer coordinate system, the inclination leads to 



(3.77) 
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and the radial systemic velocity of the observer, Vg, yields 

dx* 



dx^ 



Is 



IsVs 



(3.78) 



where 7^ = (1 — Vg)~^^'^. Although the emitter's 4-velocity has a non-zero ^-component, we 
do not need to know it, since the photon's worldline lies entirely within the 6 = tt/2 plane, 
which completely suppresses the occurance of d^/drje in the redshift formula. 

The term arising from coordinate transformation, dp/dip\e, can be evaluated, and a te- 
dious, but straightforward calculation yields 



dp 
dp 



cos cos C , 



(3.79) 



This is a standard result given by the relative orientation between the observer and the 
observed galaxy [91] and thus, is independent of relativistic effects. It is now possible to 
specify the total redshift of any emitter in an inclined galaxy with systemic velocity Vs- 
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(3.80) 
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Inserting the affine parameter dx = gtt dt gives 

/r7|:e-*{.e) + ^^^cos^cosC 



^ _|_ SVr dr 



Since the observer is far away from the origin, flat space may be assumed, so that ^ 
and the coordinate speed of hght — > 1. The denominator then simphfies to 



(3.81) 



Is (1 - Vs) = 

Finally, analogous to (3.55), the exact result is 

r 1 



l + Vs 



(3.82) 
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±\b\yhW{r) 



l + Vs 



y/1 - r^'{r) U*M 



cos^ cosC 



(3.83) 



Please note that b is only determined by (3.66), when particles on the apparent major axis are 
observed. In this case, at the moment of emission, dr^v = for the light ray, which was the 
requirement for (3.66). Particles that don't lie on the apparent major axis need a different 
mapping between b and r. 

For an exact mapping, the null trajectory of the light ray would have to be known, which 
in turn requires the metric functions <^(r) and m{r) to be known. As a first approximation, 
however, it is satisfactory to assume that the light ray's trajectory in the galactic plane is only 
bent marginally by the weak gravitational field. Using this "near flat space" approximation, 
I assume 

b = e-*(^'^--'") rN,mm ~ e-*('' '''' r sin^^r sin <f (3.84) 

as one would naturally from the classical definition of the impact parameter (fig. 3.4). It 
follows then that the observed Doppler redshift is 



Z cos ijj cos C 



Y^r <I>'(r) 
y^l - r <I>'(r) 



sin if cos ip cos C • 



(3.85) 



One of these three angles, however, is redundant and in general, only two angles are necessary 
to express the angular dependence. The relations (3.75) and (3.76) can be used to express 
the angular dependence in any of the equivalent forms 

ibsinV'cosV' / ^ . , 

sm ip cos ip cos C = —^^^^=^^^= = ± V cos^ ip — cos^ if = cos Tp sm tp . (3.86) 

Y sin^ Ip + tan^ ( 

The result is then most conveniently written in terms of the inclination angle ip and the 
angular position within the galactic disk (p, 



Z cos tp cos C 



cos ip sin (f . 



(3.87) 



y^l - r $'(r) 

With Z <^ 1 being a valid approximation, one can conclude that r<I>'(r) <^ 1 and thus 

Z cos tp cos C \/ r <I>'(r) cos tjj sin if = f Ar(r) cos ip sin ip , (3.88) 



where f7v(r) = y^r <I>'(r) was defined analogously to the rotation velocity in Newtonian grav- 
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ity. When applying the approximations z± <C 1 and Vg <^ 1 for a Taylor expansion of the left 
hand side of (3.83), it can be simplified using definitions (3.56), (3.57) and (3.88): 

1 + z± « 1 + Zg ± v^i^r) cosTp siiKf + Vs- (3.89) 

Please note that the distinction between the approaching and receding half of the disk is 
now absorbed into the position angle <f and the it in front of the Doppler term is used to 
distinguish the direction of the disk rotation field, relative to the galaxy coordinate system: 
+ for counterclockwise and — for clockwise rotation of the stars in the galaxy. Apart from 
the term Zg, this equation is identical to the equivalent in Newton dynamics [91]. 

The assumption Vg <^ 1, however, may not be applicable since there are now observations 
of rotation curves for high redshift galaxies [86], up to z ~ 1. In these cases, the correct general 
relativistic calculation must be carried out in an expanding FLRW background metric. This 
is a task that goes well beyond the scope of this thesis. 

When analyzing observational data, Z, Vg and ip are relatively easy to determine, but 
direct observations of Zg are theoretically only possible along the apparent minor axis of 
an inclined galactic disc, where s'mip = 0, i.e. no Doppler shift occurs. Practically, the 
gravitational redshift Zg is negligible compared to the Doppler redshift Z, as was shown in 
(3.70), and furthermore it is buried in noise, which generally arises from a finite thickness of 
the galactic disc. Typical observed Doppler redshifts are Z w 250^/c ^ 1 and therefore by 
(3.70), we can approximate Zg <C Z, leading to 

z± ^ ±VN{r) costp siiKp + Vg. (3.90) 

This is the same redshift relation as in Newtonian dynamics for emitters in a galaxy with 
inclination angle ^ [91]. Generalizing this result, we see that it is safe to assume that no 
general relativistic effects need be considered when determining the Doppler-shift. This is 
especially easy to understand since I have shown that the gravitational redshift Zg is always 
negligible compared to the Doppler redshift Z, and since a systemic redshift occurring from 
the expansion of the universe would affect all observed photons in the same way. 

Thus a rotation velocity curve v{r), that was obtained using common procedures in New- 
tonian dynamics, is a good approximation for the general relativistic Doppler redshift Z{r) 
and the corresponding gravitational potential <I>(r). 



3.2.5 Spacetime metric of flat rotation curve region 

As I argued in §3.2.1, the direct evidence given by observations of individual rotation curves 
suggests that the rotation velocity in the observed region is approximately constant and 
hence, implies that the dark matter density falls off as p oc r~^: 

Vc ~ const. =^ /3 oc . (3.91) 

Furthermore, observations have shown that the rotation velocity is small compared to the 
speed of light, 

250^/cw 8 X 10"^ « 1, (3.92) 
in geometrical units. From the definition 
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it follows that the energy of such a circular orbit is given by (3.39) with dr = 0: 
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In this approximation, the geodesic equation for circular motion, (3.42), yields 

,2 



d^> 

dr 



r 



(3.94) 



(3.95) 



which is identical to the corresponding formula arising in the Newtonian case. This can be 
easily integrated to 



$(r) -^Q^vl 



(3.96) 



where $o ^^nd '"o are constants of integration that have to be chosen by appropriate boundary 
conditions. This gives a good approximation to the metric component gtt for the region of 
flat rotation curves: 



gtt{r) 



=2*(r) 



-e2*o I - 



2vi 



-e2*o I - 



10" 



(3.97) 



The very small exponent 10 ^ has the effect that the metric component is almost constant 
for a very large range of values for r/rg. To see this, let —gtt e~'^'^° be bounded by 



1 -e < 
(1 - ef' < 



ro 



10- 



Then for e.g. e = 10 ^, 



10" 



-43 



< 



< 10 



,43 



< 1 + e 
<(l + e)^ 



g2*(r) _ ^2*0 



(3.98) 



Since it is unlikely that the region of flat rotation curve extends to infinity, and since there 
have been indications of a decreasing rotation curve for large radii (see §3.2.1), the integration 
constants $o 'ind rg can be determined from matching an appropriate metric at the "far-away 
radius" tq. 



Matching onto the Schwarzschild exterior metric 

The Schwarzschild exterior metric, 

ds2 = - (^1 - dt^ + (^1 _ dr^ + , (3.99) 

represents the asymptotically flat geometry that is induced by a point mass M in an oth- 
erwise empty space. Furthermore, by Birkhoff 's theorem, all asymptotically flat spherically 
symmetric vacuum spacetimes must be of this form. Thus, the Schwarzschild exterior metric 
is suitable as a far-away metric for the galaxy metric, and (assuming cosmological effects are 
negligible) can be matched with the smooth matching conditions 



S'tt, Schwarzschild (ro) 



(3.100) 
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-T' <7«,Schwarzschild(?'o) , 



(3.101) 



which yield 




2M 
1 

2M 



(3.102) 



(3.103) 



r, 



.2 ' 




SO that eliminating AI gives 




1 




1. 



(3.104) 



The result is the gu component of the galaxy metric for the flat rotation curve region, trun- 
cated at ro, merging onto an asymptotically flat exterior Schwarzschild metric: 



Similar derivations can be found in [7, 64]. 

Depending on the cosmological context, one might wish to consider different metrics as the 
exterior to the region of flat rotation curves, e.g. the de Sitter metric to represent empty space 
with a cosmological constant. Please note that the former derivation of observed redshift z± 
assumed asymptotically flat space and hence, this derivation would not hold in that form 
when using de Sitter space, which is not asymptotically flat. There might be other metrics of 
a simple form, that are suitable to model flat rotation curves and/or the asymptotic far-away 
region. 

The Tolman IV solution (see §4.3.3), for example, exhibits a density profile that goes as 
p kr~'^ + const for very large r while the pressure goes as p ~ /cr"^ — const. Thus, it 
seems that it represents flat rotation curves and asymptotic de Sitter behaviour and hence, 
might be a suitable metric for the outer region of a galaxy in a cosmology with cosmological 
constant. Unfortunately, the two free parameters of the Tolman IV solution are fixed by the 
given rotation velocity and density of the cosmological constant. The resulting distance scale 
~ 10^ kpc = 10^ Gpc is far too large to be useful for galaxies. 

3.2.6 Interpretation of rotation curve measurements in general relativity 

While the flat rotation curve implies a metric of approximate form (3.105), a rotation curve 
of arbitrary shape determines the gravitational potential <I>(r) or equivalently, the metric 
component gtti^), as was shown in §3.2.2. However, the rotation curve does not by itself 
yield any information about the metric component grr{f), i-e. the metric function m{r). 

In Newtonian gravity, the gravitational potential <I>(r) alone is sufficient to determine the 
gravitational field completely. Furthermore, it is linked to the matter density through 



which then entirely determines the matter distribution. 

In general relativity, the metric (2.32), which represents arbitrary, spherically symmetric, 
non-vacuum spacetimes and was chosen as a model for the galactic halo, is determined by 
the metric components gtt{f) and grri'''), or equivalently <I>(r) and m{r). Hence, knowledge 
of <I>(r) alone does not constrain the spacetime completely and thus, the matter distribution 




(3.105) 



V2^>(r) = 47rp(r) , 



(3.106) 
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cannot be determined uniquely. 

This can be remedied by either refining the halo model, i.e. introducing further assump- 
tions, or obtaining more data through a different observational procedure which does include 
information about grr{f)- 

Most commonly, a further assumption is made. The CDM paradigm suggests that the 
equation of state for dark matter is p{p) ~ and therefore, the source equation of weak field 
gravity (3.5), 

= 47r (p + 3p) , (3.107) 

reduces to its Newtonian counterpart (3.106). Then, the matter distribution is again given 
by the gravitational potential $(r) only. If, however, the equation of state was different, the 
derived matter distribution would of course be invalid. 

The scientifically more conservative approach is to see if there are independent means of 
measuring the equation of state. For the problem at hand, this can be achieved by observing 
galactic events that do not exclusively depend on gtt{f), but also on grrif). As I will outline 
in the next section, gravitational lensing is a suitable physical event to measure the dark 
matter equation of state in combination with rotation curve measurements. 

3.3 Gravitational lensing 

Gravitational lensing is another measurable effect of a gravitational field. Eddington first 
used the gravitational lens of the sun to test Einstein's theory of general relativity during the 
solar eclipse of the 29**^ May 1919 [27, 28]. Today, microlensing events, which temporarily 
magnify the apparent magnitude of a star, are further cases of gravitational lensing by stellar 
mass objects. 

Stronger gravitational lenses require more mass, and there are indeed observations of 
strong gravitational lensing by galaxies and galaxy clusters. Strong lensing can produce 
multiple images of the source galaxy (see §3.3.3). Lensing by galaxies on a larger distance 
scale is used to map the gravitational field of whole galaxy clusters, or using statistical 
methods, to construct a model of the gravitational field of a typical average galaxy. This 
technique is called weak lensing, see §3.3.3. 

The significant difference between measurements of the gravitational field with rotation 
curves and gravitational lensing is the speed of the probe particles. For rotation curves, 
the probe particles are the stars or the gas in the galactic disk which rotate at a clearly 
subluminal speed Vc <C c. In gravitational lensing, however, the probe particles are the 
photons of background objects, whose trajectory is bent by the foreground object. Photons 
travel at the speed of light and therefore, their perception of the gravitational field is quite 
different, compared to subluminal particles. 

In this section, I will first solve the equations of motion for a photon in a weak gravitational 
field, apply this result to gravitational lensing of point masses and then go on to a fully general 
relativistic description of lensing by arbitrary geometries using Fermat's principle and an 
effective refractive index. Finally I will comment on the actual observational techniques that 
are employed today. 

3.3.1 Solving the equations of motion directly 

The simplest case of gravitational lensing is the one where a light ray is bent by the weak 
gravitational field of a point mass. Such a lensing event is characterized by the observable 
shift in apparent position on the sky, i.e. the deflection angle Aip (see fig. 3.10). One way to 
find this deflection angle is to solve the equations of motion, as they are given by the geodesic 
equations for null curves (2.8), directly. 



46 



CHAPTER 3. GALACTIC DARK MATTER HALOS 




Figure 3.10: Deflection of a light ray by a point mass. The solid curved line demarks the 
light ray's trajectory. The angle (p{r) is the azimuthal angle for a point on the trajectory at 
radius r. the total deflection angle is given by Aip. The radius of closest approach is marked 
by r^in- The light ray's trajectory is symmetric with respect to rmin- 



The observer and the photon source are assumed to be far away from the lensing object, 
such that both can be approximated to rest in flat space at r = oo. Then, in the plane of the 
light ray's trajectory, the total change of a photon's azimuthal angle ip is twice the change 
from r = oo to r = r^i^ (fig. 3.10). If the trajectory was a straight line, this would be just vr 
and the deflection angle would be Aip = 0. Hence, the deflection angle is usually defined as 



Aip = 2 \ip{oo) - ip{r^in)\ 



vr . 



(3.108) 



The following derivation of Aip is similar to the one found in Weinberg [92, §8.5]. However, a 
new formula (3.126) is given that relates the weak gravitational field of arbitrary geometries 
specified by $(r) and m(r) to the deflection angle Aip. 

To calculate the angular difference that corresponds to the difference in radial distance 
between r = oo and r = rmin, I first combine the definition of a photon's angular momentum 
(3.32) with (3.36) which yields 







-Bir) 



~f /drV 



-4 \ d(/3"J ^ A{r) 



(3.109) 



Before integrating this equation, please note that for the radius of closest approach, r 
where dr = 0, this equation gives 



21 



A{r^ 



g2*(r-niin) 



(3.110) 



The constants of motion vanish upon inserting this into (3.109), 



dr 

d(^ 



B{r) 



or in terms of the metric functions <l>(r) and m{r), 



^('"min) 

A{r) 



(3.111) 



dr \ ' 



I 1 



2m(r) 



32'I>(r^in)-2$(r) _ 



(3.112) 
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The differentials can now easily be separated and formally integrated with the upper limit 
r = oo: 



(^(oo) — ip{r) = ± 



roc ^ 
Jr Jl - 25 



1 



df 



^1 - 2m{r)/f V(f/rmin)2 e2*{nnin)-2l>(f) _ i f 



(3.113) 



Up to this point, no approximations have been made, and the angle of deflection (3.113) is 
valid for any static spherically symmetric spacetime. Assuming that the gravitational field is 
weak in the relevant region of the galaxy, I can certainly Taylor expand the first square root: 



1 



m(r) 



m\ 2 
r 



Y^l — 2m{f)/f 

The second square root needs to be looked at more carefully. Again for weak fields: 

1 _ 

1 



(3.114) 



(3.115) 



V(f/rmin)2 - 1 + 2[cD(r^i„) - ^{f)] {f/r,,^,,Y + ' 
which can be further expanded to yield 



^(r-min) - ^{f) 
(f/rmin)^ - 1 



+ • 



+0 



(3.116) 



where the smallness of the series term still has to be shown. This term can be rearranged to 

2 



1 \r 



^(f) - ^(rmin^ 
1 - (?^min/'f)^ 



< 1 for f » Tn 



(3.117) 



However, at f = rmin both numerator and denominator vanish. Application of de I'Hopital's 
rule gives 

^'^Z. 1 - {r.^./rf = .'^Z. WU^) = ^ ^ ^'^"^ ' ' ^'-'''^ 
which has been shown to be small for the relevant region of a typical galaxy on page 41. When 
the gravitational field is reasonably smooth, both numerator and denominator of (3.117) are 
monotonically decreasing. Hence, I conclude that 



Vf > r„ 



$(f) - <I>(rn 



1 - {rmm/fY 

Finally, the second square root of (3.113), that is (3.116), can be written as 



1 



1 



1 + 



<^(f) - $(r: 



mm; ~2 



min 



$2 



(3.119) 



(3.120) 
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Therefore, the Taylor series of (3.113) for weak gravitational fields is 

1 , m{f) ^ $(f) - $(rmm) -2 

df 



2 



+0 



7] 



r 



The first term can be integrated analytically: 



1 df . ^ / mill 

arcsm 



(3.121) 



(3.122) 



The deflection angle (3.108) for gravitational lensing in a weak gravitational field is then 
given by 

Since rmin cannot be observed directly, I want to rewrite this result in terms of the impact 
parameter (3.66), 

r-min = h e*('-in) = 6 + o [$] . (3. 124) 

To the required order of accuracy, this gives 

A, . 2 /"^ ' I + ^ , (3.125) 

or equivalently 

A,.2.r^|!:i©+ia:#)4. (3.120) 



lb \/f2 — 52 

This formula is valid for all spherically symmetric spacetimes in which the lensed photons 
experience only weak gravitational fields. However, to evaluate this integral, one has to make 
a specific choice for the metric functions <^(r) and m{r), as the only quantity that is accessible 
through lensing oberservations is 

One can find an analytical solution for the Schwarzschild metric, which represents a 
"point" object with total mass M. Similarly, if the halo falls off sufficiently rapidly, the 
Schwarzschild solution will be a good approximation to the geometry outside of the halo. 
The metric components of the Schwarzschild solution (3.99) are 

m(r) = M; $(r) = ^ In ^1 - , (3.127) 

for which the integral (3.126) still cannot be solved analytically, but at least an approximate 
solution can be obtained for weak gravitational fields. Therefore, I assume that the trajectory 
of the photons passes the point mass far away from event horizon at r = 2M. Using the 
smallness of the term 

2M 2M 

< < 1 (3.128) 

I can approximate 

, . M , , 

$ r w , 3.129) 

r 

so that 

, /■°° df ( M M(r-b)] 
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Introducing the dimensionless quantity x = f/b yields the solvable integral 

2Mb dx { I 11 , , 

which evaluates to 

2M , , 4M 
A^«— {1 + 1} = — . (3.132) 

Although this method of deriving the deflection angle is straightforward, the integral 
involved can usually not be solved analytically for more general metrics. Also, the physical 
insight gained for the relation between the metric functions, <^(r) and m(r), and the deflection 
angle A(y9 is limited. To understand the influence of the metric functions on the bending of 
light of spherically symmetric objects better, it is helpful to follow the analogy of an effective 
refractive index, which relates the bending of light by gravitational fields to familiar problems 
in classical optics. 



3.3.2 Lensing with an effective refractive index 

To introduce the notion of an "effective refractive index" associated with a gravitational field, 
I will relate the path of a null curve in a static metric to Fermat's principle of least time or 
equivalently, least optical length. Then I will show how the path of a light ray in a static 
spherically symmetric spacetime is described by one scalar function only, namely the effective 
refractive index n(r). Finally, I emphasise how one can use this refractive index to easily read 
off the relation between the metric components, gu and and the deflection angle Ay?. 



Null geodesies in spacetime and the 3-dimensional Fermat metric gab 

Following Misner, Thorne & Wheeler [60, exercise 40.3], I will show that null geodesies in a 
static gravitational field can be expressed as geodesies in a Riemannian 3-space with Fermat 
metric Qah and therefore satisfy Fermat's principle of a light ray's least optical length^^. A 
static gravitational field is characterized by dtQap = and gtb = gat = 0. Hence for null 
curves with 

ds^ = g^y dX^' dX"" = (3.133) 

follows 

dt^ ={dX*f = -^dX''dXK (3.134) 
9tt 

The non-vanishing Christoffel symbols of the first kind are: 

Tttc = \d,gtt (3.135) 
Tatt = da gtt (3.136) 

Tabc = ^ [Ob gac + dc gab - Oa gbc] ■ (3.137) 
Splitting the geodesic equations (2.8) with affine parameter x into time and space yields, 

= gu^ + 2Tu.^^- (3.138) 
dx dx dx 



^^Fermat's original formulation required the integral f dt to be minimal. This statement was later shown 
to be physically incorrect and replaced with the requirement that the first variation must vanish locally, i.e. 



be extremal: 5 J dt = 0. [46] 
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d'^X'^ dX* dX* dX^ dX" 

dad _ o ^ ^att — T ^ ^abc ' 



dx^ 
d^X'i 



gad 



dx dx dx 

dX'' dX^ 



^abc — r, 



att 



9bc 
9tt 



dx dx 



(3.139) 



where (3.134) was used in the last hne. Reparametrizing with the time coordinate t = X* 
gives 



d^i dX" fdt 



0; 



2vd 



gad 



d'X 



+ 



abc att 



gtc 
gtt 



dx^ dX" 



dX'^ dh/dx' 



dt dt 



-gad ■ 



dt (dt/dx)2- 



(3.140) 
(3.141) 



Please note that the time coordinate t is not an affine parameter in the 4-dimensional space- 
time Eliminating the second derivative d'^t/dx^ from the last two equations, I find 



2 vd 



gad ■ 



d"X 



+ 



gbc 



gab 



^abc — Tatt 2 Fttc 

gtt gtt 



dX^ dX'' 



dt dt 



(3.142) 



which upon inserting the Christoffel symbols (3.135) to (3.137) takes the form 



2 vd 



gad ■ 



d'X 

dt2 



+ 



n 1 ,1 gbc r^ gab o 

Oc gab - i:(Jagbc + — Oa gtt Oc gtt 

2 2 gtt gtt 



dX^ dX' 



dt dt 



0. 



When defining the so called Fermat metric in three space dimensions, 

gab 



gab 



-gtt 



(3.143) 



(3.144) 



one can easily identify (3.143) as the geodesies of the three dimensional gab with affine pa- 
rameter t: 



2 Yd 



gad ■ 



d"x 

di2 



+ 



dc (jab ~ 2^a. gbc 



dX" dX" 



dt dt 



0. 



(3.145) 



By the basic notion of geodesies in Riemannian metric spaces (§2.2.1), this corresponds to 
minimal travelling time t in the Riemannian 3-space or analogously minimal lapsed time 
coordinate X* = t in the Lorentzian 4-spacetime: 



dj dt = 5j y^gab dX» dX* = . 



(3.146) 



This is the exact mathematical formulation of Fermat's principle, which states that the first 
variation of the lapsed time along the path between two points of a light ray's trajectory must 
vanish. An even more general result can be obtained for conformally stationary spacetimes, 



ds^ = e^/^*'^") 



{dt + (/)., dX' f + gijdX' dX^ 



(3.147) 



where exp[2f{t,X"')] is the conformal factor and <j)i{X'^) is called the Fermat one-form [71]. 
Fermat's principle of the least lapsed light travel time then takes the following form: 



Vgabdx- dXl'-4>adX" 



0, 



(3.148) 
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which reduces to (3.146) in the conformally static case with (pi = dih, where h[X'^) is a scalar 
function of the space coordinates only. 

Since the speed of light in transparent media, Vn, is given by the vacuum speed of light, 
c, divided by the refractive index n, one can formulate Fermat's principle equivalently as the 
principle of shortest optical length 

c dj = n{x, y, z) do- = , (3.149) 

where a is Euclidean arclength and 

^ = vn = -. (3.150) 
at n 

If the Fermat metric is conformally Euclidean, i.e. gab = n^(x, y, z) 6ab, equation (3.146) takes 
exactly the form (3.149), 



djn VKb dX'^ dXb = 5j n{x, y, z)da = (3.151) 



and therefore, the scalar effective refractive index n can be conceived from the metric coef- 
ficients of conformally static spacetimes g^i, that have a conformally Euclidean space-space 
part gab. 

If the space-space portion of the spacetime is not conformally Euclidean, one can still 
define an effective refractive index tensor Uab = Vg^, analogously to a refractive index 
tensor in anisotropic media like e.g. crystals [13]. This is the strong field generalization of 
the static weak field refractive index tensor presented in Appendix A. 



Non-Euclidean spaces and the propagation of light 



Although I already showed how null curves correspond to Fermat's principle, I want to present 
the relation between geodesic motion in a non-Euclidean 3-space and the propagation of 
electromagnetic waves in Euclidean space, as it is shown by Kline & Kay [46, §V.A]. 
Consider the non-Euclidean 3-space given by the differential line element 



y, z) (dx^ + dy^ + dz^) = du^ , 



(3.152) 



where n{x, y, z) is the refractive index and do" is Euclidean distance. From the laws of optics, 
we know that the propagation of light rays is governed by the relation 



d^ 



n ■ 



df\ 
d^j 



Vn. 



(3.153) 



where f = (x; y; z) is the position vector of a point on the light ray's trajectory (for reference 
see [13, §3.2.1] or [46, §V.4]). This can be reparametrized in terms of 



dr 


= n 


dr 


da 







_d 
d^ 



df\ 
"d^j 



-— n —— = n —— n —— = V n 



dx 



1 



dr \ 



and 



dxy 



(3.154) 
(3.155) 



where the gradient V is still to be taken in the flat Euclidean space. Applying the product 



52 



CHAPTER 3. GALACTIC DARK MATTER HALOS 



rule to (3.155) gives 
so that inserting 
and 

results in 



dn df d^r 
2 — — + n- 



dx ' dx^ 

dn dn dx dn dy dn dz df p 
dx dx dx dy dx dz dx dx 



1 



n^ 



df 



dx 



d^f _|_ 2 dr ^ \ df 
dx^ Vdx y dx 



df 



dx 



Vn = 0. 



Rewriting this in the usual index notation with f = Xa yields 

d^Xa , / dX'\ dXa ^ dX^ dX'' 



dx' 

which when put into the form 



dx / dx 



dx dx 



0, 



2d^Xa 



n f r 1 , \ dX^ dX^ 

n- -T— 5- + 2n < n c dab - 7; ^^bc } 3 — 

dx^ I 2 J dx dx 



0, 



(3.156) 
(3.157) 
(3.158) 
(3.159) 

(3.160) 
(3.161) 



is easily identified as the geodesic equation (2.19) of gab = ^ah with the affine connexion 

1 



Tabc — 2 {dab.c ~\~ 9ac,b 9bc,a} 



Therefore, light travels along a trajectory with 



S dx = S nda = cS dt = 



(3.162) 



(3.163) 



which has been shown in §2.2.1. Since dt = da/vn is the light's travel time, Fermat's principle 
is satisfied. 

In the next section I will show how this geometric formulation of the refractive index 
can be used to relate the metric components gu and g^r of a static spherically symmetric 
spacetime to the bending of light. 



Refractive index of a static spherically symmetric spacetime 

I will now derive an explicit formula for the exact refractive index of a spherically symmetric 
spacetime. If a conformally static spacetime is spherically symmetric, it can be expressed 
in isotropic coordinates where all information about a light ray's trajectory is contained in 
one scalar function, the effective refractive index n(r). I follow the approach of Perlick [71] 
for the exact analytic formulation and then introduce a weak field approximation that casts 
n{r) in a simple and familiar form that allows straightforward interpretation of the relation 
between the metric components and the bending of light. 

The spacetime (2.32) can always be written in conformal static form. 



ds' 



-dt^ + 



-2${r) 



1 — 2m{r) jr 



dr2 + r2e-2*M dO^ 
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= e2*W [-dt2 + 52(r)dr2 + i?2(^)d02] . (3.164) 
To perform a coordinate transformation to isotropic coordinates, 

ds^ = e2*(^-) [-dt^ + n\r) (df^ + f2 ^[^2)] , (3.165) 

one finds the conditions 

n{f)df = S{r)dr; n{r)f = R{r), (3.166) 
so that ehminating n[r) gives the transformation for the r-coordinate in integral form 

dr 



df 
f 



S{r) f 
dr ^ In — 



r a/1 — 2m{r)/r 



R{r) rc 

where rc is a constant of integration. The refractive index is then given exactly by 

R{r) 



n{r) 



f 



(3.167) 



(3.168) 



where r = r(f) still has to be determined from integration of (3.167). Although this is a simple 
and straightforward relation, the coordinate r and the refractive index can only be obtained 
analytically if the function m(r) is known. Before I go on to introduce the approximation of 
weak gravitational fields, I note that the change to isotropic coordinates allows the use of one 
refractive index n(r) only. For arbitrary coordinates, two refractive indices are necessary in 
spherical symmetry, as shown by Atkinson [2], and with less symmetry one generally needs 
an entire 3x3 tensor of refractive indices, see Appendix A. 

For weak gravitational fields, one can approximate 2m[r)/r <^ 1 and therefore, the inte- 
gral in (3.167) simplifies tremendously 



In^ 

fc 



^ 1+ 

r 



m{r) 



+ 



In h 

rc 



r 

m{r) 



2m 

r 



dr + O 



2m Y 



(3.169) 



Once one chooses the constants of integration to satisfy rc = rc, one gains the relation 



r exp 



m(r) 



dr + O 



2m 

r 



(3.170) 



Since f 



0[2m/r], the refractive index to second order in 2m /r takes the form 

2 



n[r ~ r) 



exp 



-^>(r) 



m{r) 



dr + O 



2m Y 
r J 



2m 



, , f m(r) /2?n\^ 2m o 

l-$(r) - / — ^dr + O (— J ,—$,$2 



(3.171) 
(3.172) 



The physical interpretation of this last line is simple and enlightening: 

1. The trajectory of a null curve in a static spherically symmetric spacetime is entirely 
determined by one scalar function n(r) only. This is especially also true for strong 
gravitational fields, where n(r) is given by (3.168). The requirement of two refractive 
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indices in other coordinate representations apart from isotropic coordinates is a mere 
reflection of coordinate artefacts. The physically relevant information about bending of 
light in static spherically symmetric spacetimes is encoded in a single refractive index. 

2. The refractive index can be used to calculate photon trajectories in Euclidean 3-space. 
The resulting trajectories are null curves in the isotropic coordinate representation 
(3.164) of that spacetime. Especially for asymptotically flat spacetimes this means 
that asymptotic observables, that were calculated using the refractive index n(r), have 
the same value for the Euclidean 3-space and for the spacetime. For example, the 
asymptotically observed deflection angle Aip can directly be calculated from n{r). 

3. The metric functions <I>(r) and m(r) both contribute to the refractive index and therefore 
to the bending of light. This means, that observation of gravitational lensing reflects 
contributions from both the metric components gtt and grr, unlike observations of e.g. 
rotation curves that yield information about gu only. Furthermore, if one uses the 
Newtonian approximation $'(r) « m{r)/r'^, it is evident from (3.172) that both metric 
components contribute in equal amounts to the bending of light in weak gravitational 
fields. 



It is now evident that combined observations of rotation curves and gravitational lensing 
in principle completely determine the geometry, if it can be modelled by a static spherically 
symmetric spacetime. Before I show how this can be quantified, let me illustrate how to 
calculate the deflection angle of a lensing event of a point mass and comment on current 
observational techniques. 



Derivation of Aip using the effective refractive index n{r) 

Finding the deflection angle for arbitrary spacetimes can only be achieved numerically, no 
matter which approach of calculation one persues: solving the equations of motion directly 
or utilizing the effective refractive index n{r). However, if one expects the deflection angle to 
be small, one can approximate the deflection angle analytically. 

Let e = dr/dcr be a unit length tangent vector of a light ray in Euclidean 3-space. The 
(small) total angle of deflection is then given by the difference between the asymptotic tangent 
vectors ein and Cout of the incoming and outgoing light rays, 

a = eout-ein, (3.173) 

so that for small angles a, 

Aif « |a| = ^2 Vl - cosa f« a. (3.174) 
Equation (3.173) can also be written in integral form: 

Q = J ^d(j, (3.175) 

with the integral to be evaluated along the photon's trajectory. The geodesic equation (3.161), 
or equivalently (3.153), yields 

n — = Vn — e • Vnj = V±n , (3.176) 
where Vl is deflned to give the gradient only in the directions perpendicular to e. Hence the 
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total angle of deflection is given by [77] 

d = j^^da^ jv±nda, (3.177) 

where the variation of n was approximated to be small for small deflection angles: 

ri(r) « 1 + e(r) ; with e(r) < 1 . (3.178) 

The integral in (3.177) is to be taken along the light ray's trajectory. For weak gravitational 
fields, e.g. the ones described by the refractive index (3.172), the deflection angle will be very 
small and hence, the path of integration can be approximated by a straight line, instead of 
following the exact and unknown trajectory of the light ray. 

Let the z-coordinate be pointing into the direction of this straight line of integration and 
let b = \{bx,by)\ = b\b\ be the impact parameter whose unit vector b lies in the x-y plane 
by definition. Then, both V±n and a will approximately lie in the x-y plane as well and 
by symmetry point into the b direction. Therefore, I can write for a spherically symmetric 



b 














' 



Figure 3.11: Geometry of the integration path to calculate the deflection angle. The z-axis 

forms a perpendicular triangle with the vector of the impact parameter b = bb and the 
position vector f = rr ; |6| = |f| = 1. 

refractive index n(r), 

/oo 

b-V±n{r)dz. (3.179) 
-oo 

Then, 

b ■ V±n(r) = b ■ Vn(r) =b-f^ = - n'ir) , (3.180) 

dr r 

where I used the perpendicular triangle (fig. 3.11) 

r2 = 62_^22 (3.181) 



to deduce 

From (3.181) also follows that 
so that 



b.f = cosi} = -. (3.182) 
r 



rdr = zdz (3.183) 



r b f°° n'(r) 

A(^«2 / -n'{r)dr = -2b -=AJ=dr. (3.184) 

Inserting the derivative of the refractive index (3.172), 

n'{r) ^ -^>'(r) - , (3.185) 
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gives 

This formula looks similar to the one obtained from solving the equations of motion directly 
(3.126), and indeed, the identity of both formulae can be shown by an integration by parts: 



$'(r)dr /-"^ [$(r) - $(6)]'dr $(r) - ^>(6) 



b 



oo 



Noting that the surface term vanishes after applying de I'Hoptial's rule for r = b, 1 have 
shown that the deflection angles calculated from both independent approaches, solving the 
equations of motion directly (3.126) and using the refractive index (3.186), are consistent 
with each other. 

Therefore, it is no surprise that (3.186) yields the familiar solution when inserting the 
metric functions of the Schwarzschild solution (3.99): 

^>'(r) ^ ^ ; m{r) = M , (3.188) 
into the formula for the deflection angle (3.186), 

Please note that half the deflection angle originates in the gu component of the Schwarzschild 
metric, while the other half is contributed by the grr component. This is also the case for 
the deflection angle derived via the equations of motion (3.126) and more generally for the 
effective refractive index of the Schwarzschild exterior geometry: 

M f M ^ 2M , , 

n{r)^l + /— dr = l + . (3.190) 

This "50:50" behaviour is a special feature of the Schwarzschild solution and not typical for 
more general spherically symmetric geometries. Especially when gravitational lensing occurs 
within a galaxy, where the point mass model ceases to be an appropriate approximation, 
it is necessary to realise that there are two distinct contributions to the deflection angle or 
equivalently to the refractive index. In the remainder of this section, I will outline how to 
interpret observations of lensing by extended objects. 



Advantages of the notion of the effective refractive index 

One might ask, what is the advantage of the derivation of the deflection angle using the 
refractive index? After all, it was necessary to introduce several approximations, and then 
the only result that could be obtained analytically is the deflection angle for the lensing of 
the weak gravitational field of a point mass M, which can be achieved more easily by solving 
the equations of motion directly, as in §3.3.1. 

However, while deriving this result using the refractive index, we learned how the metric 
components contribute to the lensing. Roughly speaking, for point masses, the effects are 
split up 50:50 between the gu and g^ component. We also learned that the refractive index 
n(r) is the only observable function that can be measured by observing gravitational lensing 
events. Hence, only if one assumes a-priori the quasi-Newtonian refractive index 



n(r) = 1 - 2$iens(r) 



(3.191) 
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will one measure the potential $iens(^)) which is generally interpreted as the gravitational 
potential <^(r). From (3.172), we now know that actually 

-^dr. (3.192) 

Since the only observable is n(r), it is not possible to separate the contributions from $(r) 
and m(r) by observing gravitational lensing events only. A different means of obtaining 
information about $(r) and/or m(r) is necessary, and as I have previously outlined, rotation 
curve measurements are suitable to determine <I>(r) only. 

Before I quantify the interpretation of combined observations of rotation curves and grav- 
itational lensing, I will give an overview of the current observational situation of strong and 
weak gravitational lensing. 



3.3.3 Observational techniques of gravitational lensing 

The notation in this chapter differs slightly from the previous notation to match the formulae 
as they are usually employed in the observational gravitational lensing literature. All vectors 
and gradients are two dimensional and live in the "lens plane" which will be defined later on. 



Lensing formalism as employed in the astrophysics community 

The standard formalism [10, 76, 77], which is actually used in virtually all current observa- 
tions, is based on the calculation of the deflection angle Aip using the Schwarzschild metric 
of a point mass M, as derived in (3.132) or equivalently (3.189): 

A « ^ I , (3.193) 

where a is now the total angle of defiection (previously A(p) and ^ is the 2-dimensional 
position vector of the lensed image in the lens plane. Since the distance between the source, 
lens and the observer each are all larger than the region in which most of the lensing takes 
place (~ 1^1 = 6), a thin screen approximation is generally employed. This means that all 
matter density, which is responsible for gravitational lensing, is projected into a 2-dimensional 
plane perpendicular to the line of sight. If one abandons the notion of a point mass and allows 
for a finite matter distribution, the surface density is given by the summation of all mass 
along a line of sight: 

S(6= / p{i,z)dz, (3.194) 



where S(^) is the surface density and z is the space coordinate along the line of sight. The 
defiection angle induced by a small mass dM in a small area d^^ is then given by 

da.lMi!!«|, ,3.195) 

SO that the total defiection angle of an image at position ^ is taken to be the superposition 
of all small defiections as they would occur by the isolated small masses dM: 
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This is the central formula, on which the whole formalism is based. Note that this formula 
is not restricted to spherical symmetry, but allows arbitrary matter distributions in the 
lens plane. Basically all inferences that have been made from observations of gravitational 
lensing events assume that this formula accurately describes gravitational lensing for weak 
gravitational fields. Since superposition is a valid tool in linearized weak field gravity, there 
is no objection against this approach and many consistent observations seem to support 
the validity of (3.196). However, to interpret gravitational lensing results in a first post- 
Newtonian context, one needs to be aware that the use of (3.196) implicitly assumes for a 
point mass 

r M 

^'lcns(r) = $(r) = / ^dr, (3.197) 

which can be concluded from comparison of (3.190), (3.191) and (3.192). When applying the 
previously introduced superposition, it is straightforward to see that the spherically symmet- 
ric Poisson equation for Newtonian gravity implies 

A$iens(r) = 47rpiens(r) fiens(r) = J ^^^^^ dr . (3.198) 

As I argued previously in §3.2.6, this is equivalent to assuming that the pressure of the matter 
which generates the gravitational field vanishes or is negligible (|p| <C p)- If this assumption 
is not imposed a-priori, gravitational lensing oberservations will only yield information about 
$iens(^)- The actual mass distribution m(r) can then only be obtained from (3.192) when 
^>(r) is known from a different kind of observation, e.g. rotation curves. Hence, miens('') and 
Plensii") are the spherically symmetric matter and density distributions under the assumption 
that pressure does not contribute significantly to gravitational lensing. 

Without derivation, I note the three basic equations of gravitational lensing, following 
from (3.196), that relate the source and image properties [10, 47, 76]: 

1. Arrival time: 

T = \ (ii - - i^iii) . (3.199) 

The appropriately scaled arrival time r of an image at position S^j is defined by this 
equation. ^5 is the position of the source of the image. The time delay between two 
images is given by the difference of their arrival time r. The lensing potential 
is the appropriately scaled 2-dimensional equivalent of the 3-d gravitational lensing 
potential: 

^P{iI) oc J «5iens(e/, ^) dz a 2 y S(|) In \ij - || d^^ (3.200) 

2. Lens equation: 

6 = 6 + W(e/) • (3.201) 

The lens equation is equivalent to Fermat's principle of least arrival time, Vr = 0, and 
gives the position of the source for one of its images at ^/ if the lensing potential ^p 
is known. 



Magnification tensor: 



, .-1 d^T _ __d^^d^ 



K + fj. 
K — fi 



(3.202) 



This equation defines the inverse of the magnification tensor My at the image position 
The components of the magnification matrix are given in suitably rotated coordi- 
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nates of the lens plane which diagonalize the matrix. The total magnification is given 
by M = ^detMj"^^ = (k^ ~ /"^) ^- denotes the expansion of the image and fx 
denotes its shear. 

The arrival time r and the lensing potential ip are suitably scaled to include all relevant 
cosmological information about the distances and expansion between the observer, lens plane 
and the source [76]. The lensing potential is defined by the lens equation to satisfy 

Vij = a. (3.203) 

Therefore, a spherically symmetric mass distribution yields by (3.196) [47, 77] 

9^^(^) = l!!!i£££M, (3.204) 
r 

where m\cns{r) is the inferred mass within a sphere of radius r under the assumption that 
pressure in the matter which causes the bending of light is insignificant. 

Before I give a brief overview over the use of these formulae in observational data analysis, 
I want to point out that these three equations depend on the lensing potential tp or its first 
and second derivatives, all of which are related directly to <I>icns and its first and second 
derivatives. 



Strong Lensing 

The term strong lensing is generally used to describe gravitational lensing that happens 
within 1-10 kpc of the lensing object, which usually leads to the occurrence of multiple images. 
Despite the name "strong" lensing, the gravitational fields involved are weak enough to allow 
the use of the linearized formalism that was introduced earlier. At present there about 70 
multiple image systems of galaxies or galaxy clusters known^^ [49], with more data being 
collected and therefore increasing this number. Due to the small number of images, typically 
2-4, the number of observable parameters - image positions, amplification and time delay 
- is also very limited. Furthermore, conclusions about the lensing potential ip{£,i) are only 
possible at the few image positions. Hence, lensing systems that have their multiple images 
at different radii are especially suitable for probing the potential ^*ions(^) or the mass profile 

"^lens('^)- [47] 

The data from a strong lensing system is mostly fitted with some promising parametric 
model for the mass distribution, such as e.g. the singular isothermal sphere (which corre- 
sponds to a flat rotation curve for all radii), or the NFW halo [62]. This procedure allows 
a reasonable mass estimate, despite the limited number of data points, but it assumes the 
shape of the mass profile miens ('") a-priori [47]. 

Recently, with the upcoming availability of high-resolution observations of lensing sys- 
tems, non-parametric methods are being developed and employed [25, 26, 74, 93]. These 
methods allow the reconstruction of a "pixelated" density map of the surface density S from 
the lens equation (3.201) and hence, are suitable for actually measuring the mass profile 
iTT'icnsi'i") without any prior assumptions about its shape. 

While the time delay (3.199) between the lensed images was initially thought to be a 
suitable tool for measuring the Hubble parameter Hq [47, 50], it turned out that the time 
delays are very sensitive to the mass distribution [65] and estimated values for the Hubble 
parameter tend to be too low [50] when compared to cosmological measurements of the Hubble 
parameter (see §3.1.5). Conversely however, if an accepted value for the Hubble parameter 

^''For an up to date list see http://www.cfa.harvard.edu/glensdata/. 
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is assumed, time delay measurments yield important constraints on the mass distribution. 
Unfortunately, time delay mesurements are still very rare [50], since extensive monitoring of 
the image flux is necessary to determine the time delay between the images which ranges 
from a few days to years. 

Weak Lensing 

Weak lensing involves even weaker gravitational fields than strong lensing. Generally the 
images of source galaxies are lensed by the weak gravitational field of a foreground galaxy, at 
a distance of ~ 10-1000 kpc. Because the distance between lens and image is so much bigger 
than for strong lensing, the number of source/lens galaxy pairs for weak lensing is very large. 
On the downside this also means that there are no multiple images or significant distortion 
and hence, the signal of a single weak lensing event is practically not detectable. However, 
due to the abundant occurence of weak lensing events one can constrain the halo potential 
of an "average galaxy" by combining a statistically large sample [15]. 

Since there are no multiple images, time delay is not accessible for weak lensing and the 
position of the image is also not measurably different from the source position. Hence, the 
only way to obtain a signal from weak lensing is by analyzing the small distortion of the 
images that turns spherical sources into ellipses that are tangentially aligned with respect to 
the vector that connects the lensed image with the center of the lensing potential well. A 
large number of lensed images then should draw up a coherent pattern of the gravitational 
potential in a galaxy field [15]. Of course the majority of lensed galaxies are not spherical in 
their unperturbed appearance, but already elliptical. The complex image ellipticity e, often 
also called the "image polarization", contains the magnitude of the ellipticity of an image 
and its orientation. For a completely random distribution of unperturbed images of source 
galaxies, one would expect the ensemble average to vanish: (Sunpcr) = 0. The weak lensing 
signal hides in the fact that the observed average ellipticity does not vanish: (eobs) 0- 
Therefore, the extraction of the weak lensing signal from a large sample involves rigorous 
statictics which are explained in detail in e.g. [34, 41, 45, 44]. 

Straightfoward averaging over all lensed images produces a picture of an "average galaxy" 
that is not necessarily characteristic of any individual galaxy. To understand the gravitational 
potential of galaxies in their diversity, one has to apply appropriate scaling laws in the weak 
lensing data analysis [47]. The first critical scaling is the distance of the lens and the source 
galaxies. This means that ideally the data contains the redshift of all involved galaxies. If 
redshift is not available, distance statistics have to be applied which reduces the quality of 
the data [15]. The second scaling is to assume a universal density profile that is followed by 
all galaxies. This is again a parametric approach to determine the mass profile miens (j")! as 
was already discussed for strong lensing [47]. Hence, weak lensing does not yield an unbiased 
mass profile that is free of prior assumptions. 

The last issue of weak lensing that I want to mention in this context is that weak lensing 
can only use (3.202) to determine the second derivative of the lensing potential ip. Therefore, 
only the first derivative of the mass distribution can be inferred. The constant of integration 
that is necessary to obtain an absolute value for the mass profile miens (j^) is not always 
unambiguous, and hence the mass profile can generally not be uniquely determined [47]. 
Recent efforts to break this so called mass-sheet degeneracy combine both weak and strong 
lensing techniques [14]. 

3.3.4 Interpretation of gravitational lensing observations 

As I have outlined in §3.3.2, gravitational lensing observations of a spherically symmetric 
weak gravitational field described by (2.32) are entirely determined by the scalar refractive 
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index 

n(r) = 1 - 2«>iens(r) . (3.205) 

Different types of lensing observations yield information about <l>icns('') and/or its first and 
second derivative. $iens('") is usually interpreted as the Newtonian gravitational potential 
which leads to the conclusion that the density and mass of a spherically symmetric object 
are given by 

P\cns{r) = A$icns(?') ; micns(?^) = r'^ ' (3.206) 

However, I already showed that if one does not make this assumption about <I'icns('^)i but 
instead interprets it in terms of the metric functions <l>(r) and m{r) for spherically symmetric 
weak gravitational fields, one finds the relation (3.192), 

^dr, (3.207) 

or equivalently, 

o d<I> 

2 miens (r) ~ + '"(^) • (3-208) 

This relation, which handles non-negligible pressure in a static spherically symmetric weak 
gravitational field correctly, is the key equation for combining measurements of rotation 
curves with gravitational lensing. 

3.4 Combined measurements of rotation curves and gravita- 
tional lensing 

The simplest metric that can be used as a realistic model for a galaxy is a static spherically 
symmetric spacetime as it is given by (2.32): 

1 — 2m{r)/r 

Naturally, this metric stops to be an appropriate description if the considered galaxy is too 
aspherical or not dominated by a spherical dark matter halo. For all other cases, the system 
is completely described by the two metric functions $(r) and m(r). 

As I have shown in §3.2 and §3.3, neither rotation curve nor gravitational lensing obser- 
vations on their own can determine both functions uniquely. Moreover, in both cases it is 
necessary to make an assumption about the equation of state, p = p{p), of the matter which 
generates the gravitational field. Without such an assumption, the data cannot be inter- 
preted. The commonplace assumption is p = for which the interpretation of both kinds of 
observations takes the familiar Newtonian form of weak gravitational fields. 

Along the way, I indicated that combined observations of rotation curves and gravitational 
lensing eliminate the need for such an a priori assumption and furthermore, allow to deduct 
the equation of state - or depending on the quality of the data, at least limits on it. I will now 
quantify how one can interpret measurements of combined rotation curve and gravitational 
lensing observations. 

3.4.1 Formalism combining rotation curves and lensing 



The relation between the gravitational potentials that are obtained by rotation curve or 
lensing observations, <^rc and ^icns) and the metric functions $(r) and m(r) has been derived 
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previously in §3.2.4 and (3.192): 

^Rc(0 = $(r), (3.210) 
$lcn.(r) = -$(r) + - y -^dr. (3.211) 

In general, these potentials do not agree. Only in the Newtonian limit, where the pressure is 
negligible, are they identical: <I>rc = ^icns = *^*N- Since this is the standard assumption for 
interpreting rotation curve and lensing data, most observations are processed to report the 
mass profile instead of the gravitational potential - which also avoids the problem of finding 
a suitable constant of integration for the potential. Therefore, it is useful to discuss the 
combined measurements of rotation curves and gravitational lensing in terms of the inferred 
mass distributions directly. For both kinds of measurements, the usual assumption is that 
the Newtonian relation 

mncnUr) - ^^^^^^ = <I>;,,/,_(0 (3.212) 

defines the mass distribution in terms of the obtained potential. Under this assumption, the 
relations between the obtained "pseudo-mass" distributions and the metric functions are: 

mRcir) = r2$'(r) (3.213) 
miens(r) = ^r^<^'{r) + ^m{r). (3.214) 

These "masses" describe the observations equivalently to the potentials (3.210) and (3.211). 
The quotation marks around "masses" indicate that these profiles no not generally describe 
the total mass-energy within a sphere of radius r. Recalling (3.5), 

V2$Rc(r) = V2$(r) ^47r{p + pr + 2pt) , (3.215) 

and the M-Einstein equation (2.38), 

m(r) 



47rp(r), (3.216) 



one can easily deduce the Poisson equations 

V2$Rc(r) « 4TT{p + pr + 2pt) (3.217) 

V2$iens(r) ^ 4tt (^p+^\pr + 2pt]^ . (3.218) 

The physical significance of the measureable "mass" distributions rriRc ^-nd miens becomes 
now clear by integrating these Poisson equations: 

mRcir) ^ i-rr [p{r) + pr{r) + 2pt{f)] f"^ df (3.219) 
Jo 



pr 

?Tilens('^) ^ 4:7r 

Jo 



p{r) + ]^Pr{f) +Pt{f) 



f^df. (3.220) 



It is now evident that the Newtonian assumption (3.212) only leads to its intended physical 
interpretation as total mass-energy within a sphere of radius r when the sum of the principal 
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pressures vanishes: 



m(r) m^c{r) ^ micns{r) 



Pr{r) + 2pt{r) w . 



(3.221) 



Therefore, I shah refer to these "masses" , as they are determined by the estabhshed rotation 
curve and lensing formahsms, as pseudo-masses. Inverting (3.213) and (3.214) determines 
the metric functions in terms of these observed pseudo-masses: 



$'(r) 
m(r) 



2 miens (?^) - mRc{r) 



(3.222) 
(3.223) 



The gravitational potential <^(r) can be obtained by integration under suitable boundary con- 
ditions. However, only the derivatives of <I>(r) contribute to the field equations and hence to 
the density and pressure profiles. The exact Einstein equations (2.38)-(2.40) for a spherically 
symmetric stress-energy tensor of the form (2.66) are: 



Stt p(r) 
8'irpr{r) 
87rpt(r) 



2m' (r) 



m{r) 



r^>'(r) ^1 



2 m(r) 



-— ^m' (r) r — m{r)^ [l + r$'(r)] + 



+ 1 



2 m{r) 



$'(r) 



+ $'(r)2 + $"(r) 



(3.224) 
(3.225) 

(3.226) 



Inserting (3.222) and (3.223) and their derivatives gives the density and pressure profiles of 
the galaxy in terms of the observed pseudo-mass functions mRc('") and miens (''): 



47r r'^p(r) 
Att r'^pr{r) 

Att r'^pt{r) 



„ rriRcir) - micnsir) 



+ 



2 m 
r 



mRc(r) - miens (?-) 



[A-Kr'^Pr{r)]' + 



- 1 




+ 









2 m 

r 



(3.227) 
(3.228) 

(3.229) 

(3.230) 



As consistency checks, I note that: 

• The Einstein equations in curvature coordinates, (3.227)-(3.229), agree to the given 
order of 2m/r with the Einstein equations of the metric in isotropic coordinates (3.165), 
and therefore the approximation f ~ r as of (3.170) is valid; 

• From the equations (3.227)-(3.229) follows that 

A7:r^[p{r)+Pr{r)+2pt{r)\=m'^c(T)+0 [ — ] , (3.231) 

and 



47r 



m 



lens 



'2m 

2 m 

r 



(3.232) 
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and thus these results are consistent with the weak field approximation of the field 
equations (3.5), and the interpretations of the pseudo-masses (3.219) and (3.220); 

Should the observed mass profiles agree with each other, i.e. m = ttirc = "^lonsi the 
density and pressure profiles yield the desired result of the Newtonian limit: 



Airr"^ p{r) 



m! 


ir), 


O 




^2m 






O 


/2m 






y r 



(3.233) 
(3.234) 

(3.235) 



I conclude that the currently existing formalisms for analysing data from rotation curve 
and gravitational lensing observations can be used to obtain the pseudo-masses m^c{r) and 
'^icns(?')) which by (3.227)-(3.229) yield the density and pressure profiles in a first post- 
Newtonian approximation. Furthermore, from the combination 

47T r2 (p, + 2pt) « 2 {m'^c " "^Us) , (3-236) 



one can immediately infer that the observed system is Newtonian, in the sense of negligible 
pressure content, if and only if m'^Q{r) « "T-iens('")- Furthermore, defining the dimensionless 
quantity 

^ Pr{r) + 2pt{r) _ 2 m^Rc(^) - mlnsi^) 

Mr) ^ 3 2m[^^^ir)-m'^^{r) ^" 

gives a convenient parameter that determines a "measure" of the equation of state. 



3.4.2 Observational situation 

The post-Newtonian formalism I have outlined requires the simultaneous measurement of 
(pseudo-) density profiles from rotation curve and gravitational lensing observations. 

While in principle these profiles do not have to be of the same galaxy, they must be 
comparable in the sense that they accurately describe "similar" galaxies. For example, weak 
lensing measurements (§3.3.3) can be used to statistically infer the (pseudo-)density profile of 
an "average" galaxy [15]. At the same time, analysing the dynamics of satellite galaxies gives 
the rotation curve and thus, the corresponding pseudo-density profile, of another "average" 
galaxy [15]. Whether these two "average" galaxies are comparable or not depends on many 
factors, such as e.g. the distribution of galaxy morphologies in both samples, the statistical 
noise, the employed models for the (pseudo-)density distribution, etc. These statistical issues 
render the fast-growing collection of weak lensing data problematic for applying the formalism 
presented. 

On the other hand, combined simultaneous measurements of rotation curves and lensing 
of individual galaxies are extremely well suited for our formalism. However, while there 
is a large number of individual rotation curves available (> 100,000; [86]), the number of 
individual "strong" lensing systems with multiple images is rather limited (~ 70; [49]^^). 
Combined observations are further aggravated by the differing distance scales: Most high 
quality rotation curves are naturally available for galaxies with a low to intermediate distance, 
corresponding to a redshift of up to z ~ 0.4 [86], while gravitational lenses are easier to detect 
at intermediate to high redshifts of z > 0.4 [49], since the image separation scales increasingly 

^*For an up to date list see http://www.cfa.harvard.edu/glensdata/. 
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with the redshift of the lensing galaxy [76, 50]. Therefore, even for nearby galaxies with 
existing combined measurements of kinematics and lensing (e.g. 2237+0305 at z ~ 0.039 
and ESO 325-G004 at z w 0.035), the lensing data is resticted to the core region, while the 
rotation curve is only described by few data points in the outer region of the lens galaxy 
[4, 85, 88]. Consequently, the inferred pseudo-mass profiles are available for the same galaxy, 
but unfortunately at different radii and therefore not comparable. 

Although the observational situation makes it currently difficult to employ the formalism 
presented, the situation is likely to improve in the future when observations with a higher 
resolution will be carried out ~ preferably with an emphasis on obtaining high-resolution 
rotation curves for lensing galaxies that exhibit lensed images at different radii. 

For instance, detailed observation of the recently discovered closest known strong lensing 
galaxy ESO 325-G004 [85] could provide satisfactory data to allow the decomposition of 
density and pressure of the galactic fluid, as outlined in this article. The system consists 
of an isolated lensing galaxy at redshift z ~ 0.035 with an effective radius of Res = 12"5 
and arc-shaped images of the background object at ~ 3", and possible arc candidates at 
R ~ 9". Smith et al. intend to collect more detailed data [85] that hopefully will include 
extended stellar dymanics and hence, allow for a direct comparison of the rotation curve and 
lensing data, if the arc candidates at i? ~ 9" turn out to contribute to the measurements. 

3.4.3 Conclusions 

In this chapter, I have argued that the standard formalism of rotation curve measurements 
and gravitational lensing make an a priori Newtonian assumption that is based on the CDM 
paradigm. I introduced a post-Newtonian formalism that does not rely on such an assump- 
tion, and furthermore allows one to deduce the density- and pressure-profiles in a general 
relativistic framework. In this framework, rotation curve measurements provide a pseudo- 
mass profile mRc(^) and gravitational lensing observations yield a different pseudo-mass 
profile w,iens(^)- Combining both pseudo-masses allows one to draw conclusions about the 
density- and pressure profiles^^ in the lensing galaxy. 



In the case of absent or negligible pressure, this could be used to observationally confirm the 
CDM paradigm of a pressureless galactic fluid. Conversely, if significant pressure is detected, 
a decomposition of the galaxy morphology would allow new insight into the equation of state 
of dark matter. 

Since the formalism presented is based on a first-order weak field approximation, I suggest 
that to confirm the findings, one should re- insert the obtained density and pressure profiles 
into the metric (3.209). The actual observed quantities can then be extracted numerically 
for comparison from the exact field equations (3.224)-(3.226) and the geodesic equations. 

Even though data might not yet be available to constrain the dark matter equation of 
state noticably, one should note that the possibility of non-negligible pressure in the galactic 
fluid introduces a new free parameter into the analysis of combined rotation curve and lensing 
observations. Therefore, the approach presented here might also help to shed some light on 
prevailing problems that arise when combining rotation curve and lensing observations. For 
example, an unresolved issue exists when measuring the Hubble constant from the time delay 
between gravitationally lensed images: Using the standard models for matter distribution in 

^^These formulae are given in SI units, hence the factor of (? . 




(3.238) 



(3.239) 
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the lens galaxy, the resulting Hubble constant is either too low compared to its value from 
other observations, or the dark matter halo must be excluded from the galaxy model to obtain 
the commonly accepted value of Hq [50]. A possible explanation of this trend might lie in a 
disregarded pressure component of the dark matter halo. 

Finally, I want to emphasize that the key point of this new formalism is that in gen- 
eral relativity, density and pressure both contribute to generating the gravitational field. 
Furthermore, the perception of this gravitational field depends on the velocity of probe par- 
ticles. These effects become especially important when one compares rotation curve and 
gravitational lensing measurements, where the probe particles are fundamentally different: 
interstellar gas or stars at subluminal velocities for rotation curves, and photons which travel 
at the speed of light for lensing measurements. The presented formalism accounts for these 
crucial differences between the probe particles, and relates observations of both kinds to the 
the density and pressure profile of the host galaxy. Although only static spherically symmet- 
ric galaxies are considered in this thesis in a first post-Newtonian approximation, the basic 
concept is fundamental and can be extended to more general systems with less symmetry. 
A suitable framework for considering most exotic weak gravity scenarios is provided by the 
effective refractive index tensor, as introduced by Boonserm et al. [12], see appendix A. 



Chapter 4 

Compact spherical objects 



After using the spherically symmetric spacetime (2.32) to investigate general relativistic ef- 
fects in galaxies, I now want to focus on compact objects with an extremely high density. 
While the necessity of general relativistic effects for galaxies depends on the pressure content 
in the galactic fluid, it is certain that the strong gravitational fields of compact objects cannot 
be described adequately by Newtonian gravity. 

4.1 Compact stars and black holes 

Today, basically three types of compact objects that require general relativity for their de- 
scription are well established in the astrophysics community [60, §24]: 

White dwarfs are remnants of not so massive stars that used up their nuclear fuel. They 
no longer have a source of energy and radiate away the remaining thermal energy. The 
pressure gradient that holds up a white dwarf arises from the degeneracy pressure of 
electrons. Their maximum mass is given by the Chandrasekhar limit of ~ 1.4 A^© [18]. 
General relativity does not have a significant influence on the overall stellar structure, 
but is necessary to describe the pulsation frequencies accurately. 

Neutron stars are remnants of more massive stars that are held up by the degeneracy 
pressure of neutrons. Therefore they have approximately the density of atomic nuclei. 
General relativity is necessary for the appropriate description of the stellar structure, 
due to the strong gravitational fields inside the neutron star. 

Black holes are the only possible configuration in general relativity for an object composed 
of ordinary matter, once it becomes more compact than the Buchdahl-Bondi bound [56] 
of 2M/R < 8/9. They exhibit a curvature sigularity at the center and an event horizon 
at the Schwarzschild radius of Rs = 2M. Rotating and charged black holes have a 
few more peculiar properties that shall not be discussed in this work. Black holes do 
not exist in Newtonian gravity and hence, general relativity is indispensible for their 



While white dwarfs and neutron stars only require general relativistic corrections, I now 
want to concentrate on black holes which are only conceivable in general relativity. Static 
black holes (without electric charge) are generally described by the Schwarzschild metric 



description. 




dr^ 



(4.1) 



1 - 2M/r 



an 
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which was the first solution of Einstein's field equations (2.33) to be discovered by Karl 
Schwarzschild in 1916 [81]. Its mathematical simplicity is also its physical curse: The solution 
exhibits a curvature singularity at the center where r = 0, which is considered to be a 
physically pathologic feature that hints that general relativity alone might not be enough 
to describe the physics of such strong gravitational fields in such a small region. Maybe a 
unified theory of gravity and quantum mechanics will provide the necessary solution. Until 
such a solution is found, physicists are happy to accept that the pathological singularity is 
"hidden" behind the event horizon at Rs = 2M, which acts as a "one-way membrane" from 
which not even light can escape, and anything that passes this horizon is lost in the black 
hole without hope of recovery^ . 

Other than the central singularity, which cannot be detected in principle, the event horizon 
is the only directly detectable feature of a black hole. The indirect observational evidence 
for black holes - through the gravitational fields surrounding it - has recently become very 
convincing, as I already outlined in §3.1.2. However, until there is direct evidence of the event 
horizon, the indirect observation only tells us that there is a very compact object, and the 
only widely accepted object that fits this category is a black hole. 

In this chapter, I will explore basic properties of mathematical alternatives to black holes. 
The use of the word "mathematical" indicates that the presented models satisfy all math- 
ematical requirements of forming a compact object with a radius of i? ~ Rg^ while the 
physical reality of the underlying assumptions is currently as observationally unaccessible as 
the physical implications of the inside of a black hole. 



4.2 Gravastars 

One of the small number of serious challenges to the usual concept of black holes is the 
"gravastar" (gravitational vacuum star) model that was recently developed by Mazur and 
Mottola [58]. In the gravastar picture, or the very closely related quantum phase transi- 
tion picture developed by Laughlin et al. [19, 20], the quantum vacuum undergoes a phase 
transition at or near Rg where the event horizon would have been expected to form. 

In the Mazur-Mottola model, a suitable segment of de Sitter space^ is chosen for the 
interior of the compact object while the outer region of the gravastar consists of a (relatively 
thin) finite-thickness shell of stiff matter {p = p) that is in turn surrounded by Schwarzschild 
vacuum (p = p = 0). Apart from these three explicitly mentioned layers, the Mazur-Mottola 
model requires two additional infinitesimally-thin shells with surface densities a± , and surface 
tensions 'd±, that compensate the discontinuities in the pressure profile and stabilize this 5- 
layer onion-like construction, effectively by introducing delta-function anisotropic pressures 
[58]. Since infinitesimally thin shells are a mathematical abstraction, for physical reasons it 
is useful to minimize the use of thin shells, for example by attempting to replace the thin 
shells completely with a continuous layer of finite thickness, like in [17]. 

This section is a brief summary of the results that were found by Cattoen, Visser and 
myself [17]. The full article can be reviewed in appendix B. 



^Attempts to combine black hole physics and quantum theory have led to models in which black holes emit 
"Hawking radiation" [9, 36] and maybe even information that was previously thought to be lost forever [37]. 

^De Sitter space has an equation of state p = —p > 0, which is also identified as the equation of state of 
the cosmological constant and hence, the interior of a Mazur-Mottola gravastar is sometimes said to be filled 
with "dark energy". 
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4.2.1 Generalized model without thin shells 

In the article [17] we demonstrate that the pressure anisotropy imphcit in the Mazur-Mottola 
infinitesimally thin sheh is not an accident, but instead a necessity for ah gravastar-hke 
objects. That is, attempting to build a gravastar completely out of perfect fluid will always 
fail. (Either the gravastar swells up to infinite size, or a horizon will form despite one's best 
efforts, or worse a naked singularity will manifest itself.) We derive this result by working 
with configurations where pressure is assumed continuous and differentiable, and analyzing 
the resulting static geometry using first the isotropic TOV equation (2.77), and then the 
anisotropic TOV equation (2.76). 

To have a useful model, we should retain as much of standard physics as possible, while 
introducing a minimum of "new physics". In the spirit of Mazur and Mottola [58], and 
Laughlin et al. [19, 20], we will keep the density positive throughout the configuration but 
permit the pressure to become negative in the gravastar interior. To avoid infinitesimally thin 
shells one must then demand that the radial pressure Pr is continuous (though the density 
need not be continuous, and typically is not continuous at the surface of the gravastar). 
Qualitatively the radial pressure is taken to be that of fig. 4.1. 



Pc 



core crust atmosphere 














/ 






^ 

r 













Figure 4.1: Qualitative sketch of a gravastar, explicitly labelling the "core", "crust", and 
"atmosphere", refers to the radius at which the local gravitational acceleration g = 
^ r (p + 3pr)/(l — 2m/ r) vanishes; for details see §B. 

That is, to "smooth out" the infinitesimally thin shells of the Mazur-Mottola gravastar 
model, we shall consider static spherically symmetric geometries such that: 

• Inside the gravastar, r < R, the density is everywhere positive and finite. 

• The central pressure is negative, Pc < 0, and in fact Pc = —Pc- 
(We do not demand p = —pr = —pt except at the centre.) 

• The spacetime is assumed to not possess an event horizon. 
This implies that Vr we have 2m(r) < r. 

These three features, positive density, negative central pressure, and the absence of horizons, 
are the three most important features characterizing a gravastar. They also ensure that a 
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gravastar does not exhibit the pathological central sinuglarity and the event horizon of a 
black hole. Instead, the spacetime metric of a gravastar is well-defined everywhere. 



4.2.2 Bounds on the pressure anisotropy 

It is straightforward to demonstrate that a gravastar-like object with negative central pressure 
cannot consist entirely of a perfect fluid. First, I rewrite the anisotropic TOV equation (2.76) 
as 

dpr _ 4:TT r {p + Pr) {p + 3pr) ^ 2pA 

dr 3 1 — 2m{r) jr r ' 

where I made use of the average density p{r) = m(r)/(|7rr'^). The anisotropy parameter A 
is defined as 

A ^ , (4.3) 

P 

and hence, for A = 0, (4.2) represents the isotropic TOV equation. 

To realize the failure of a perfect fluid in a gravastar, note that the maximum pressure in 
an ordinary spherical compact object occurs at r = 0, where dpr/dr = naturally for A = 0. 
In a gravastar, this extremum is a pressure minimum and the maximum pressure occurs at 
a certain radius r^ax 7^ 0, see fig. 4.1. Therefore, at rmaxj the first term of the right hand 
side of (4.2) is negative for a gravastar, while the left hand side vanishes by the definition of 
a pressure maximum. Consequently, A must be positive, i.e. anisotropic pressure is required 
to allow a pressure maximum for r > 0. A more detailed discussion of the region in which 
perfect fluids fail to deliver a valid solution can be found in appendix B. 

Once we accept that perfect fluid spheres are not what we are looking for to model 
gravastars, one might wonder what happens to the Buchdahl-Bondi bound for isotropic fluid 
spheres. It has been shown that for p' < and pt < Pr the 8/9 bound still holds. However 
if the transverse stress is allowed to exceed the radial stress, pt > Pr-, then the upper limit 
shifts to 2M/R < k < 1, where k depends on the magnitude of the maximal stress anisotropy 
[33]. In the gravastar picture, I already showed that A > and hence, pt > Pr- Therefore, 
the compactness of a gravastar is not limited by the Buchdahl-Bondi bound, but only by 
the magnitude of the maximal pressure anisotropy and the regularity of the metric, i.e. 
2m{r)/r < 1. Let us now make these qualitative considerations more quantitative. To do 
this, let us rewrite (4.2) in the form 



^ ^ Pt-Pr ^ ^ 
P 2 



Pr_^f-^_^ Pr \ Trn + Anpr r"^ 



P \ p / [1 — 2m/r] 



(4.4) 



In appendix B we have determined the smallest interval in radii for which anisotropic pressure 
is necessary to be (r'g,rniax]- We can now ask for explicit bounds on A in that interval. 

For the interval r G [ro,rmax\, by making use of p', > and Pr > 0, we find the simple 
lower bound 

1 2m/r , , 

A > - > 0. 4.5 

4 1 - 2m/r 

Now in the region [ro,rmax] we have pt > Pr > 0, consequently if the dominant energy 
condition [DEC; p > and \pi\ < p] is to be satisfied we must at the very least have A < 1. 
But this is guaranteed to be violated whenever 2m/r > 4/5. 

That is: If the gravastar is sufficiently close to forming a horizon, in the sense that 
2m/r > 4/5 somewhere in the range [ro,rmax]) then the DEC must also be violated at this 
point. ^ Consequently, any gravastar that is sufficiently close to forming a horizon will violate 



■'And even if we were to discard the entire positive-pressure region {ro,R), we can nevertheless still apply 
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the DEC in its "crust" . 

For the interval r G (rg,ro) we find the considerably weaker bounds 

- 2p 2p il- 2m/r ^ ' 

where we have used j?', > 0, pr- < 0, and the null energy condition [NEC; p + Pi > 



4.2.3 Minimizing the anisotropic region 

Let us now attempt to minimize the region over which anisotropy is present. It is easy to see 
that at rg we have 

A(r,) = ^>0. (4.7) 
z p 

Therefore, in the case where the anisotropy is confined to the smallest interval possible, we 
want Pr{rg) = 0, corresponding to an inflexion point for the radial pressure. At the point tq 
of zero radial pressure, the anisotropy cannot vanish: 

A(ro) = ^ + i '"^Z^" >0. (4.8) 
^ 2p 4 1 -2m/ro ^ ' 

At the point of maximal radial pressure, the anisotropy also has to be non-zero, at least if 
we take the limit from below: 

Mr^j = \ . f 1 + ^Ui + ^^"i > ° . (4.0) 



4 1 - 2m/ 

^max \ P J \ ^ 

Beyond the peak Vr > rmax, it is possible to arrange A = 0, though at a price: If we wish 
to confine the anisotropy to the smallest interval possible, we have to set A(r '"max) — 
which leads to a discontinuity in p'^ and A as well as a "kink" in the pressure profile pr at 
'"max- (However, p^ itself is still continuous, as is the density p.) Indeed we then have 

pUr-iax) = A(r^,J. (4.10) 

'max 

The implications of confining the pressure anisotropy to the smallest interval possible are 
shown in fig. 4.2, where the anisotropy is confined to the region r G (rg,rmax]- 

The Mazur-Mottola model is now recovered as the limiting case where ^ tq <— Tmaxj 
and so all the important anisotropy is confined in their inner thin shell. The anisotropy 
A — > oo, because p' — > oo. Effectively A is replaced by choosing an appropriate finite (in 
this case negative) surface tension •& and surface energy density a which is given by the 
Israel-Lanczos-Sen junction conditions [5, 43, 53, 82, 90]. 

The second, outer thin shell which is present in the Mazur-Mottola model is not a physical 
necessity, but is merely a convenient way to avoid an infinitely diffuse atmosphere that would 
arise otherwise from the equation of state p = p. A finite surface radius R can be modelled by 
altering the equation of state slightly to include a finite surface density ps which is reached 
for vanishing pressure: p{p) = ps +p. Then, the outer thin shell can be omitted when joining 
the gravastar metric onto the Schwarzschild exterior metric. 



this bound at ro itself: If 2m(ro)/ro > 4/5 then the DEC is violated at ro. 
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Figure 4.2: Qualitative sketch of radial pressure as a function of r for a gravastar that is as 
near as possible to a perfect fluid. Note the inflexion point at rg and the kink at rmax- 

4.2.4 Gravastars vs. black holes 

The article [17] discusses the results outlined above in more detail. We explored the gener- 
alization of the Mazur-Mottola gravastar and found that any object that binds all its mass 
within a radius of i? ~ Rs = 2M with 2M / R > 8/9 must have anisotropic pressure in 
certain regions, i.e. pt > Pr- Therefore, the Buchdahl-Bondi bound does not apply, and the 
compactness is only bound by the absence of an event horizon, i.e. 2m(r)/r < 1. 

Consequently, a black hole and a gravastar cannot be distinguished by their gravitational 
effects only. Only a direct measurement of the gravastar surface or equivalently, the black 
hole's event horizon, would provide evidence for one model or the other. While the event 
horizon is in principle detectable through the effects of being a "one-way membrane", i.e. 
the absence of bremsstrahlung from infalling particles, a prediction for measureable surface 
effects of a gravastar is not available as yet. 

The gravastar model exists in a mathematical context, but there is no generally accepted 
physical explanation for a spherical object that exhibits the characteristic pressure profile of 
a gravastar. 

4.3 Anti-Gravastars 

Following our findings that gravastars with negative pressure near the center must have a 
crust that consists of some sort of anisotopic fluid [17], I now investigate another possible 
configuration that departs from the conventional assumption that the density and pressure 
in a stellar object are both positive. I call an object of finite extent with positive pressure, 
p > 0, everywhere and negative mass-energy density at the center, pc < 0, an anti- gravastar. 

Such an object will, of course, have a negative total mass M if the energy density p 
is negative throughout the body of the anti- gravastar. Should objects with a net negative 
mass exist, they would be observable through an unusual microlensing-profile [21] or a self- 
accelerating "-|-M/— M-rocket" [11] - both of which have not yet been observed, and probably 
never will be. Since there is no observational evidence that such negative mass objects exist. 
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we consider this type of anti-gravastar only in a mathematical fashion. 

There is, however, another type of anti-gravastar with a positive total mass M > 0. It 
consists of a negative density core with a mass of Mcore = —M- < and a positive density 
crust with Mcrust = M- + M > (see fig. 4.4). When matched onto a Schwarzschild exterior 
metric, these anti-gravastars could in theory mimic the gravitational field of a usual spherical 
stellar object. 

Unlike the gravastar, we shall see that an anti-gravastar can consist entirely of a perfect 
fluid. To see that, I assume the anti-gravastar to consist of a perfect fluid and then find that 
the isotropic TOV equation satisfies the conditions that constitute an anti-gravastar. 

4.3.1 Anti-gravastar model 

The most important properties of a static, spherically symmetric anti-gravastar are: 

• Inside the anti-gravastar, r < R, the pressure is positive and finite everywhere. 

• The central density is negative, pc < 0. 

• The induced spacetime does not have an event horizon, i.e. Vr : 2m{r) < r. 

To close the set of differential equations (4.11) that determine a static spherically sym- 
metric spacetime for a perfect fluid, 

m' = Airpr'^] (4-11) 

/ {p + p){m + A-Kpr^) 

^ r'^[l — 2m{r)/r] ' 

two initial conditions and a third equation that fixes the relation between p and p are neces- 
sary. For the initial conditions, I choose naturally 

1. ?n(0) = since I want the anti-gravastar to have a regular origin and 

2. p{R) = which defines the finite radius of the surface R where I want to match a 
Schwarzschild exterior metric with M = m{R). Equivalently, one can choose p(0) = 
Pc > and define the surface radius R by the condition p{R) = 0. If no surface can be 
defined this way, one needs to adjust the assumptions made. 

In the following sections, I will present two analytical geometries that represent a perfect 
fluid anti-gravastar and therefore prove their possible (mathematical) existence. These are 
the constant density Schwarzschild interior geometry and the less known Tolman IV solution 
[24]. Both geometries constitute anti-gravastars with a negative total mass. 

Although I wasn't able to find an analytical example of an anti-gravastar with a positive 
total mass, I will present the basic features of an equation of state that should lead to a 
perfect fluid anti-gravastar with a positive total mass. Such an object could in principle 
mimic the gravitational field of a black hole, just like the gravastar presented earlier. 

4.3.2 Schwarzschild interior anti-gravastar 

The Schwarzschild interior solution represents an object whose density is constant throughout 
its entire body. Due to its simple formulation, it is often used to illustrate basic physical 
concepts, although a constant density object is generally not thought to be an accurate 
description of astrophysical bodies. 
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In the context of anti-gravastars, I highlight the properties of the Schwarzschild interior 
solution, if one assumes negative constant density. The metric for the Schwarzschild interior 
solution is given by [24] 



A- B\ 1- ] dt'+{l-—r\ dr^ + rMQ 



C2 



C2 



(4.12) 



Calculating the ti-component of the Einstein tensor, this leads to 

3 



SttGnC^ 



> 0. 



Hence, to consider negative density I substitute the constant 
force negative density p^:. The resulting metric is: 



(4.13) 

-C^ which obviously will 



ds' 



A-B\ 1 



Q2 



dt'+{l + 



C2 



dr^ + r^dO^ 



with a negative constant density 



p* 



< 0. 



SttGnC^ 

The rr-component of the Einstein tensor yields the well-known density profile 



pir) 



1/0*1 



y^l - 2M/R - v^l - 2m{r)/r 
3^/1 - 2M/R - Vl - 2m(r)/r 



(4.14) 



(4.15) 



(4.16) 



where I found that A = 2,B^l - 2M/R, B is arbitrary and C = ±R^^-R/2M. 

Since all masses in (4.16) are negative, all square-roots exist. From m{r)/r oc r^, I 



conclude that 



yr <R: 



2m{r) 




2M 


r 


< 


R 



(4.17) 



Therefore, both numerator and denominator of (4.16) are positive. The negative density p^, 
appears in (4.16) with an accompanying minus-sign which has been replaced by the absolute 
value for clarity. Thus, it is ensured that the pressure profile is positive right up to the surface 
R. Hence, the modified Schwarzschild interior metric (4.14) represents an anti-gravastar. 

One should note that a brief calculation shows that 



yr <R: y/l^^2MjR. > a/I - 2m(r)/r 
Thus, the NEC is always violated. 



p* +p{r) < 0. 



(4.18) 



4.3.3 Properties of the Tolman IV solution 

The Tolman IV solution is another analytical geometry representing a spherically symmetric 
perfect fluid object. In its standard form, it is given by the metric [24] 

Ag2 _ ( ^ , ^ ^^2 , l + 2r2/^2 

- ^ ^^ + ^2^ + (l-r2/C2)(l + r2/^2)^^ +r au . (4.iyj 
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The pressure and density profiles follow from the Einstein equations: 

= 8.G,CHA^ + 2r^) ^ ^''^'^ 
_ ?,A^{A^ + C^) + {lA^ + 2C2)r2 + Gr^ 
^^"■^ " 87rG,vC2(^2 + 2r2)2 ' ^^'^^^ 

If we assume the constants and to be positive, we immediately see that Vr : /? > 0. 
From (4.20) we find that, if we define the surface by p{R) = 0, the surface radius is given by 
R? = |(C^ — A^) and thus, if R? > 0, both pressure and density are positive for r < R. The 
parameter B is irrelevant for the whole discussion - it can easily be absorbed into a modified 
time-coordinate - and hence, will be neglected from now on. 

The expressions for the pressure and density simplify tremendously at the center for r = 0, 

_ C^-A^ _ 3{A^ + C^) 

and therefore it is possible to express the constants and in terms of the central pressure 
Pc and the central density pc- 

^42 ^ 3 1 ^ (g2 ^ 3 1 

47r Gat pc + 3pc ' 47r Gat pc - 3pc' 

It is generally assumed that the parameters and are positive - hence the square- 
notation - but in general, both parameters can have negative values. To keep the notation 
sensible, this is realised by inverting the accompanying sign of the parameter in the metric 
and all other derived formulae. Therefore, I distinguish four cases: 

1. ^2 > 0, G^ > ^ pc> 2,\pc\ > 

2. A2<0,G2>0 <^ 3pc<-\pc\<0 

3. ^2 > 0, G^ < <^ 3pc > \pc\ > 

4. ^2 ^ 0, G^ < <^ pc< -3\pc\ < 

In the following, I will discuss the properties of the Tolman IV geometry in all four cases and 
then show which one is suitable as an anti-gravastar. 



Naked singularities. Probably the most obvious conclusion from (4.20) and (4.21) is that 
when A^ < there will be a pole in the pressure and density at = \A^ /2\ which implies a 
naked curvature sigularity (see e.g. §B) and therefore geometries of that kind should only be 
utilized in regions where r'^ ^ \ A^ /2\. In short, the Tolman IV solution with ^4^ < cannot 
cover the entire spacetime without invoking a naked curvature singularity. 

Existence of surface. Since the surface (if it exists) is located at R^ = ^{C"^ — A^), we can 
conclude immediately that there cannot be a well defined surface (in the sense of p{R) = 0) 
for case 3 and that there always will be a surface for case 2. Note that the latter will also 
have a naked singularity either below, above or even at the surface. 

To see what happens to the other cases, it is helpful to rewrite the surface radius in terms 
of the central pressure and density: 



,2 ^ 3 Pc 

2TTGNpl-{^Pcf 



(4.24) 



76 



CHAPTER 4. COMPACT SPHERICAL OBJECTS 



Thus, cases 1 and 4 will only have a well-defined surface ifi' pc > 0. 



Equation of state. Due to the simpicity of the pressure and density profiles (4.20) and 
(4.21), which are simple rational functions, it is possible to eliminate and gain the equation 
of state that is inherent to the Tolman IV solution: 

, , _ 8^^ {AttGnpC^Y + {13 A' + 2C') {AttGnpC^) +^{2A' + C^f 

AitGnCHA^ + 2C^) ■ ^^^^ 

This equation of state is quadratic in the pressure p and takes a relatively simple form when 
expressed in terms of the central pressure and density: 

( ^ ^ 2 , 5pc-llPc , {Pc-Pc){pc-3pc) , . 

P\P) = P H PH (4.26) 

Pc +Pc Pc+ Pc Pc + Pc 

{IQp^ + 5p, - llp.y , ^ 
32{pc + Pc) 

Asymptotic behaviour. The asymptotic behaviour of the Tolman IV solution is obtained 
by Taylor-expanding (4.20) and (4.21) about r — > oo: 



3 A^ + 2C^ 1 

Knlr^oo - i67rGjvC2 ^ 32ttGnC^^^ 

3 ^2 ^2^2 1 
p(r)\ = \ ■ h O 

PK )\r^oo IGttGnC^ 32TTGNC^r^ 

That is, both pressure and density approach a constant value. 



(4.27) 
(4.28) 



= = 16nLc^ ' ^'-''^ 

and both profiles exhibit a r~2-falloff close to infinity which implies that 

m[r oo) oc r + "cosmological background oc r^" . (4.30) 

In a universe with a cosmological constant, this could in principle represent a dark matter 
halo of a galaxy with a flat rotation curve. However, the density of the cosmological constant, 
the rotation velocity of halo objects and the size of a typical galaxy give three constraints 
that would have to fit the two free parameters and of the Tolman IV geometry: 

• The density of the cosmological constant can be worked out from the values given in 
table 3.1: poo ~ 0.77 x lO-^^kg/m^. Thus, it follows from (4.29) that 

(7^ PS 1.2 X 10^^ kpc^ w (350 Gpc)^ . (4.31) 

• The mass distribution in a fiat rotation curve region is typically 



m(r) 2 f 250km/s 



2 



v'rot ~ ( " ) ~ 7 X lO"'^ const. , (4.32) 



r \ c 

so that the density in the region of fiat rotation curves is given by 



P=^.^l. (4.33) 

47rr2 47r r2 
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Then, by comparison with (4.27), I conclude that 

42 _|_ 9/^2 

7 X 10"^ A^ ^ -2C^ w 2.4 X lO^^kpc^ , (4.34) 



8C2 

meaning that if the Tolman IV geometry is to describe a galaxy with flat rotation 
curves exhibiting a rotation velocity of v^-ot ~ 250km/s, embedded in a universe with 
a cosmological constant that has a density of poo ~ 0.77 x 10"'^'^ kg/m^, the parameter 
will be negative and the parameter will be positive. This corresponds to case 2. 

• Since < 0, the Tolman IV solution exhibits a naked singularity at = \A'^/2\. 
This means that the asymptotic behaviour, as it is illustrated by the Taylor expansion 
in (4.27), is only valid for r 3> |A|/\/2. That is, the Tolman IV geometry with an 
asymptotic density of poo ~ 0.77x 10"^'^ kg/m'^ exhibits the flat rotation curve behaviour 
(with the given velocity) only for 

r > 3.5 X lO^kpc, (4.35) 
which is of course a distance scale that is far too large for typical galaxies. 

Summarizing, the Tolman IV solution matches asymptotically onto either de Sitter space (for 
> 0, cases 1 and 2) or anti-de Sitter space (for < 0, cases 3 and 4). Therefore, it 
could prove to be a useful geometry if one wishes to embed a spherically symmetric object 
into a cosmology with a cosmological constant. However, the two free parameters constrain 
the geometry severely, e.g. I have shown that the Tolman IV geometry cannot represent the 
flat rotation curve region of a typical galaxy. 

As a last remark about the asymptotic behaviour, I note that the density and pressure of 
Tolman IV at infinity are related to their central values through C^: 

Poo = -Poo = ^ • (4.36) 

Derivatives of the profiles. The derivatives, with respect to r, of the pressure and density 
profiles are 

^ ~ "AttGnC^ (A2 + 2r2)2 ' ^^'^^^ 

^ {A^ + 2C'){5A' + 2r')r 
^ 47rG7vC2(^^ + 2r2)3 ' ^ ' 

From that we see that p' never changes sign and is zero only for r = and r — > co. Whether p' 
is increasing or decreasing is determined by the constants and C^: signp' = —sign [6*2(^2 + 
2C2)]. 

For ^42 > 0, the sign of p' is the same as that of p'. If yl2 < (cases 2 and 4), the sign of 
p' first changes at the curvature singularity at r"^ = —A'^/2 and then there is an extremum 
at r2 = — 5/2 ^42 where the sign also changes. Thus sign p' = signp' in the central region and 
near infinity, and sign p' = — signp' for —A'^/2 < < —5/2^2. 

Total mass. The density profile can easily be integrated to yield the mass m(r) that is 
enclosed in a sphere with radius r (Gat = 1 for convenience): 

, , (^2 + C72 + r2) r3 
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Thus, the mass grows as for r ^ cxd as one would expect from a "cosmological background 
density". Subtracting this asymptotic mass gives 

which grows as r for r — > co as I already mentioned in the paragraph about asymptotic 
behaviour: 



3 



A^ + 2C72 



(4.41) 



Although this behaviour is typical for the halo of flat rotation curve galaxies, I have already 
shown that the Tolman IV geometry cannot be parametrized suitably to match the densities 
and distance scales of a realistic galaxy. 

If the configuration has a well defined surface where p{R) = 0, the total mass M = m{R) 
is given by 

VWc^Z^_^ 

and thus we see that the total mass will be negative if is negative and if there is a well 
defined surface radius. 



Anti-Gravastar configurations of the Tolman IV metric 

We have seen that the Tolman IV solution can represent quite different geometries, depending 
on the choice of the sign of the parameters A^ and C^. It can exhibit nasty features like a 
naked curvature singularity or have a well-defined surface with vanishing pressure. In any 
case, the asymptotic behaviour of the Tolman IV solution is either de Sitter or anti-de Sitter 
space. 

A quite surprising finding was that the Tolman IV solution has a simple, analytical 
equation of state. It represents a two-parameter (^4, C) sub-class of fluid spheres that have 
a quadratic three-parameter (a, /3, 7) equation of state of the form 

p{p) = a{A,C)p^ + C)p + -i{A,C). (4.43) 

The three parameters a, /3 and 7 can be expressed through either and or pc and pc- 
Thus, if one is considering a quadratic equation of state for a spherical object, one might 
want to check whether it falls into the Tolman IV category. 

The usual object that will be considered is likely to have all parameters positive, which 
is case 1 with pc > 0. This will result in a "normal" compact object with a radius R and a 
total mass M = R-^/SC^. 

If we are looking for configurations that qualify as an anti-gravastar, we are immediately 
restricted to cases 3 and 4 with Pc > and pc < 0. Since case 3 does not have a surface, we 
can rule it out already. For case 4 we know r'^ < R'^ < /2 and thus, the naked singularity 
will occur outside the surface and therefore shall not be of concern. The derivatives of both, 
pressure and density, are negative for r < R and thus, both will be decreasing, i.e. the density 
will always be negative. Thus, the total mass M will be negative, which also can be deduced 
from < 0. 

Finally, we find that the modified Tolman IV solution for a negative mass anti-gravastar 
is ^ 

- B yi ^2 J ^ [l + r^/C^){l-r^/A^f 
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Figure 4.3: Pressure profile (left) and density profile (right) in dimensionless units for A = 5, 
C = 2 and thus R = V7 {Gn = 1, c = 1). Note that p{R) = - gg^^ is finite. 



From the negative density at the surface, we can immediately conclude that the NEC will 
be violated there. In fact, it turns out that 

(A^ + 2C^)(r^ - A^) 
since r'^ < R^ < A^ /2 < A^. Thus, the NEC is violated everywhere. 



4.3.4 Equation of state for an anti-gravastar with positive total mass 

Since we already have seen that perfect fluid anti-gravastars with negative total mass can 
exist, one might ask whether perfect fluid anti-gravastars that have a positive total mass can 
exist in theory as well. I will show that such objects can exist when one chooses the equation 
of state carefully. 

An anti-gravastar can only have a positive total mass if there is "more" positive density 
matter than negative density matter. "More" in this context means the absolute value of the 
positive mass is greater than the absolute value of the negative mass, |M_|_| > 

The simplest configuration that can have a positive total mass is a negative density core 
with a positive density crust (fig. 4.4). This satisfies our requirement that the central density 
must be negative. The positive density matter in the crust shifts the problem of converting 
infalling "normal" matter into negative density matter from the surface to inside the anti- 
gravastar. This conversion process still has to be justified in a physically relevant model, but 
for the time being, we just explore the mathematical implications. From that point of view 
it is more convenient to model the density-transition with a continuous equation of state and 
match a positive mass Schwarzschild exterior metric at the surface R. 

A simple equation of state that can in principle represent all necessary features is 

pip) = P tanh ^ + PS (4.46) 

which is plotted in fig. 4.5. The parameter po > marks the pressure for which the density 
changes its sign, D > is a measure for the width of the transition region and ps > denotes 
the surface density, which has to be non-zero if one wants to avoid an infinitely far streching 
atmosphere. 
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Figure 4.4: Pressure and density for an anti-gravastar with positive total mass. 

It is easy to check that this choice of equation of state satisfies the NEC: 

p + p = p ^1 + tanh^5^^ +P5 > P5 > 0. (4.47) 
For Pc ^ po and pc ^ I also find 

Pc = Pc tanh ^"^^"^ + P5 ~ -Pc + PS ^ -Pc (4.48) 

and thus, this equation of state corresponds approximately to anti-de Sitter space near the 
center. 




Figure 4.5: Equation of state (4.46) with po = 2, D = 0.1 and ps = 0.1. 
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To clarify the following discussion, I use the average density p and the local gravitational 
acceleration g as they were defined in §4.2. Additionally, I note that 

sign p' = -sign [{p + p) g] = -sign [{p + p){p + 3p)] . (4.49) 

The simplest anti-gravastar will have p' < throughout the interior. Then, for pc > 0, 
the pressure is guaranteed to be positive inside the surface radius R, which is defined by 
p{R) = 0. Hence, the first property of an anti-gravastar in §4.3.1 is satisfied. If pc > po, 
the equation of state (4.46) makes sure that the central density is negative, which was the 
second property of an anti-gravastar. The last important property of an anti-gravastar is the 
absence of event horizons, i.e. 2m[r)/r < 1. This property can generally only be evaluated 
numerically. Thus, to construct an anti-gravastar with the equation of state (4.46), it is 
sufficient to show that p' < everywhere when pc > Po- However, the absence of horizons 
still needs to be verified numerically. 

By (4.49) and knowing that the NEC is satisfied, I conclude that 

p <0 <^ g(xp + 3p>0 <^ p{r)>-\p{r), (4.50) 

i.e. the local gravitational acceleration needs to point inwards everywhere. Since the equation 
of state (4.46) induces a pressure and density profile with the qualitative shape of fig. 4.4, it 
is safe to assume that 

Pc < p{r) (4.51) 
everywhere. Furthermore, from (4.48), one can easily conclude that 

-Pc = \Pc\ < Pc 
for Pc > Pq. Combining these last three inequalities yields 

P{r) >^Pc>-^Pc>-^ p{r) => g>0. (4.53) 

This inequality guarantees positive g everywhere. However, it is too restrictive in the outer 
regions of an anti-gravastar. But since p > everywhere, it is sufficient to apply the less 
strict bound 

p>0 g>0 (4.54) 

for the outer regions of the anti-gravastar. I conclude that if the parameters of (4.46) are 
chosen carefully enough to 

1. ensure that pc < 0, i.e. pc > Po] 

2. have gravitational attraction everywhere, i.e. g > which is generally established by 
the inequality (4.50) and specifically by either (4.53) or (4.54); 

3. avoid the formation of an event horizon for all r < R, i.e. 2m{r)/r < 1, 

the equation of state (4.46) will lead to a perfect fluid gravastar. This equation of state can 
potentially yield an anti-gravastar with a positive total mass. To ensure this, the parameters 
of (4.46) must be chosen to also 

4. make sure that the positive density crust dominates the negative density core, i.e. 

I ^crust I ^ I ^core \ • 

Finding an actual set of suitable parameters, and therefore proving that anti-gravastars 
with a positive total mass are mathematically possible, requires one to solve the TOV equation 



4.52 
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numerically. This task was too time consuming to be realized within this thesis and therefore, 
I leave it to be carried out in the future. 

One remarkable fact that arises from this configuration, is that the density in the crust 
can be much higher than the density in an object with the same total compactness 2Ad/R. 
Since the local compactness 2m{r)/r gets a "negative headstart" from the core, the crust 
actually has to contain more positive density matter than another object with the same total 
mass M and radius R. 

4.3.5 Summary 

I have found two analytical solutions that can represent anti-gravastars with a negative total 
mass: the Schwarzschild interior solution and the Tolman IV geometry. Both solutions consist 
of a spherically symmetric perfect fluid. Among the known spherical perfect fluid solutions 
[24], there may be more that actually can represent an anti-gravastar. However, a systematic 
search has not yet been carried out. 

Unfortunately I was not able to spot an analytical solution that constitutes an anti- 
gravastar with a positive total mass, although there is no indication that such a solution 
should not exist. To emphasise the possibility of a positive mass anti-gravastar, I used the 
equation of state 

p{p) = P tanh ^ + PS (4.55) 

to illustrate the necessary conditions for the existence of such an object. However, the 
ultimate proof for positive mass anti-gravastars, using numerical methods, is yet to be carried 
out. 

While the definition of an anti-gravastar already violates the WEC and DEC near the 
center where p < 0, there also seems to be a correlation between the NEC and the sign of 
the total mass of an anti-gravastar. Both negative total mass solutions, the Schwarzschild 
interior and Tolman IV, violated the NEC everywhere whereas the example for a possible 
positive mass anti-gravastar satisfied the NEC everywhere. 

The version of the positive mass theorem obtained by Penrose, Sorkin & Woolgar [70] 
states that if the NEC is satisfied everywhere, the total mass must be positive. Applied to 
the example in §4.3.4, this theorem implies that if an anti-gravastars with an equation of 
state (4.55) exists, it must have a positive total mass. 

Conversely, one can conclude that an anti-gravastar with a negative total mass must 
violate the NEC at least locally in some region. However, it is an open question whether a 
local (or also global) violation of the NEC is permittable for a positive total mass or whether 
violation of the NEC immediately leads to a negative total mass. 

Analogously to gravastars, if anti-gravastars with a positive total mass and a sufficient 
compactness of 2M / R ~ 1 existed physically, they would exhibit an external gravitational 
field that is indistinguishable from that of a black hole. Although the model of an anti- 
gravastar is not favoured by me personally as an alternative for black holes, it demonstrates 
- just like the gravastar - that the current observational evidence for compact objects in the 
universe does not automatically lead to the conclusion that it must be a black hole — there 
is still a little wriggle room left. 



Chapter 5 

Summary and conclusions 



The discussions in this thesis have shown that pressure plays a special role in general relativity 
as opposed to Newtonian gravity. In general relativity pressure is a source of the gravitational 
field, and the perception of that gravitational field by passive particles depends on the particle 
speed (§3). Since gravity and fluid dynamics are distinct theories in the Newtonian picture, 
these "relativistic" aspects of pressure cannot be adequately described in a quasi-Newtonian 
approximation. 

5.1 Rotation curves, gravitational lensing and the equation of 
state of the galactic fluid 

After giving an overview of typical galactic structure and possible dark matter candidates in 
§3.1, I went on to discuss how to interpret rotation curve measurements in general relativity 



It has been shown previously by Lake [52], (and also by Nucamendi, Salgado and Sudarsky 
[64]), that the observable rotation curve data of edge-on galaxies contain contributions from 
gravitational redshift, as well as the Doppler shift which arises from motion of the tracer 
particle towards or away from the observer (§3.2.2). However, applying this result in an ob- 
servational context, I found that the gravitational redshift (3.70), as opposed to the Doppler 
shift (3.71), contributes at a higher order of the small rotation velocity to the observable 
frequency shift. Therefore, the gravitational redshift is negligible when processing the obser- 
vational data. Furthermore, I extended the derivation of the general relativistically derived 
total redshift to be applicable to arbitrarily inclined galaxies, and found that under the per- 
mittable weak field approximations, the resulting total redshift (3.90) is identical to that 
found in Newtonian gravity. Consequently, the established data reduction procedures for 
rotation curve measurements yield the general relativistically correct redshift for weak gravi- 
tational fields. However, as I argued in §3.2.6, the interpretation of that observed redshift in 
general relativity is different from the widely adapted interpretation in Newtonian gravity: 
The potential $rc('^) = ^(r), which can be determined from the redshift, is generally not 
equal to the Newtonian gravitational potential ^^{r), given by the field equation (3.106), 



in §3.2. 



V^$Ar(r) = 47rp{r) . 



(5.1) 



Instead, the appropriate weak field equation in general relativity is^ (3.5) 



V2$Rc(r) =4^ [p{r)+3p{r)]. 



(5.2) 



To keep the conclusion as simple as possible, I use p in the sense of p = (pr + 2pt)/3. 



QQ 
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The Newtonian limit is recovered iff \p\ <C p. 

In §3.3 I pointed out that the currently employed formalism for gravitational lensing is 
based on a superposition of the deflection angle of a point mass. Since this superposition 
excludes the notion of pressure, it cannot possibly model gravitational lensing appropriately 
when the lensing galaxy is made of some form of matter with a non-negligible pressure 
content. To find a description that allows for possible pressure in the active gravitational 
matter, I argued that the effective refractive index is the only possible observable quantity of 
gravitational lensing by a spherically symmetric galaxy. In isotropic coordinates (3.172), the 
effective refractive index is a scalar, while other coordinate systems require a refractive index 
tensor with at least two independent components (see §A). This is still true for aspherical 
configurations as long as the 3-dimensional space part of the metric is conformally Euclidean. 
By comparison of the refractive index (3.172), 

n(r) ^ 1 - 2cDi,,3(r) « 1 - ^r) - j'l^dr, (5.3) 

with its commonly used quasi-Newtonian counterpart, 

n{r) = 1 - 2$iens(r-) ~ 1 - $7v(r) , (5.4) 

I conclude that the lensing potential '^icns('^)) as obtained from the established data analysis 
methodology, has to be interpreted in terms of the metric functions <I>(r) and m(r) in the 
following manner: 

cl>iens(r) = ^$(r) + ^ l^dr. (5.5) 
Therefore, the corresponding field equation for the lensing potential (3.218) is 



V2$lcns(r-) = 47r 



3 

p{r) + -p{r) 



(5.6) 



Again, the Newtonian limit is recovered iff \p\ <C p. 

In the case of non-negligible pressure, (5.2) and (5.6) can easily be inverted to yield the 
density p and pressure p of the galactic fluid when sufficient data is available from independent 
rotation curve and lensing observations of the same galaxy. However, the critical word in 
the previous sentence is "sufficient", as it is not easy to simultaneously obtain quality data 
for rotation curves and gravitational lensing for the same object. A more detailed discussion 
of the problems involved is given in §3.4.2. The basic problem is that good rotation curve 
data at all radii can be obtained for nearby galaxies, but for these galaxies the lensing data 
tends to be concentrated towards the center. Since the recent past has shown us that the 
observational situation is only getting better, I am confident that in the near future there will 
be enough quality data available to employ the formalism presented and discussed in §3.4.3. 

It will then be possible to infer the equation of state of the galactic fiuid directly from 
observations, which might confirm the Cold Dark Matter paradigm observationally or, if 
significant pressure is detected, even shed new light on the nature of dark matter. 



5.2 Black holes vs. gravastars: Mathematical simplicity or un- 
resolved physics? 

Black holes are widely accepted in the astrophysics community because they are a mathemat- 
ically simple solution in the strong gravity regime of compact objects, that also fits the latest 
observational evidence for a very dense matter concentration at the center of the Milky Way 
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and other galaxies. However, the simpHcity comes at the price of a pathological curvature 
singularity at the center of a black hole, which indicates that general relativity indeed fails 
to deliver an appropriate description of physics at this singular point. A unified theory of 
gravity and quantum theory might resolve this issue, but it is not necessary to appeal to 
the "holy grail" of the great unified theory to find a solution that provides us with compact 
objects similar to black holes but without mathematical pathologies. 

The gravastar picture, as developed by Mazur and Mottola [58], provides a mathematical 
alternative to black holes: A massive core of matter with the equation of state of the cosmo- 
logical constant, p = —p > 0, is covered in a thin but finite thickness shell of stiff matter with 
an equation of state p = p. Implicit in this model are two infinitesimally thin shells at the 
boundaries between the core and the stiff matter shell as well as between that shell and the 
surrounding vacuum. The spacetime of a gravastar is free of curvature singularities and is 
regular in the sense that there is no horizon. Together with Celine Cattoen and Matt Visser, I 
have generalized this layered gravastar to a model that exhibits a continuous pressure profile 
without distinguished layers. We showed that gravastars without any thin shells must have 
anisotropic pressure in their "crust". Therefore, the inner thin shell is unavoidable if one 
wants to model a gravastar with perfect fluid only. 

Since the gravastar core is made out of matter similar to the cosmological constant, it is 
also said to be filled with "dark energy" , a form of mass-energy that is widely accepted in 
cosmology where the dark energy density is very low. In a gravastar, this dark energy would 
have a density that is much higher than predicted by other standard astrophysical models. 

A different, even more exotic alternative to black holes is the anti-gravastar. Matt Visser 
and I explored the properties of a compact object that has positive pressure everywhere but a 
core consisting of negative density matter. I argued that objects of this kind could also mimic 
the external gravitational field of a black hole. However, even though the anti-gravastar can 
consist entirely of a perfect fiuid, it is more exotic in the sense that it consists of matter with 
a negative density. Unlike dark energy, there are no observational indications that matter of 
this kind might exist in the universe. 

While being agnostic about the existence of gravastars, anti-gravastars or black holes, 
I want to point out that there is as yet no direct evidence for black holes with an event 
horizon and a central curvature singularity. The observations of the most compact objects 
are consistent with a black hole, but also with a gravastar or an anti-gravastar. However, 
before gravstars or anti-gravastars can be considered a serious alternative to black holes, their 
physical properties should be explored in more detail. Being a mathematical solution does 
not make the (anti-) gravastar a physically relevant object. 

I hope the research presented in this thesis will make a valuable contribution to prevailing 
questions of astronomy, and I would like to thank again all people who helped me creating 
this work. 
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Effective refractive index tensor for 
weak-field gravity 

Petarpa Boonserm, Celine Cattoen, Tristan Faber, 
Matt Visser, and Silke Weinfurtner 
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Victoria University of Wellington, 
P.O.Box 600, Wellington, New Zealand 

Gravitational lensing in a weak but otherwise arbitrary gravitational field can be described in terms 
of a 3 X 3 tensor, the "effective refractive index". If the sources generating the gravitational field 
all have small internal fluxes, stresses, and pressures, then this tensor is automatically isotropic and 
the "effective refractive index" is simply a scalar that can be determined in terms of a classic result 
involving the Newtonian gravitational potential. In contrast if anisotropic stresses are ever important 
then the gravitational field acts similarly to an anisotropic crystal. We derive simple formulae for the 
refractive index tensor, and indicate some situations in which this will be important. 
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APPENDIX A. WEAK-FIELD REFRACTIVE INDEX TENSOR 



Appendix B 

Gravastars must have anisotropic 
pressures 

Celine Cattoen, Tristan Faber, and Matt Visser 

School of Mathematics, Statistics, and Computer Science, 
Victoria University of Wellington, 
P.O.Box 600, Wellington, New Zealand 

One of the very small number of serious alternatives to the usual concept of an astrophysical black hole 
is the "gravastar" model developed by Mazur and Mottola; and a related phase-transition model due 
to Laughlin et al. We consider a generalized class of similar models that exhibit continuous pressure 
— without the presence of infinitesimally thin shells. By considering the usual TOV equation for 
static solutions with negative central pressure, we find that gravastars cannot be perfect fluids — 
anisotropic pressures in the "crust" of a gravastar-like object are unavoidable. The anisotropic TOV 
equation can then be used to bound the pressure anisotropy. The transverse stresses that support a 
gravastar permit a higher compactness than is given by the Buchdahl-Bondi bound for perfect fluid 
stars. Finally we comment on the qualitative features of the equation of state that gravastar material 
must have if it is to do the desired job of preventing horizon formation. 
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APPENDIX B. GRAVASTARS & ANISOTROPIC PRESSURES 



Appendix C 



Combining rotation curves and 
gravitational lensing: How to 
measure the equation of state of 
dark matter in the galactic fluid 

Tristan Faber and Matt Visser 

School of Mathematics, Statistics, and Computer Science, 
Victoria University of Wellington, 
P.O.Box 600, Wellington, New Zealand 

We argue that combined observations of galaxy rotation curves and gravitational lensing not only 
allow the deduction of a galaxy's mass profile, but also yield information about the pressure in the 
galactic fluid. We quantify this statement by enhancing the standard formalism for rotation curve 
and lensing measurements to a first post-Newtonian approximation. This enhanced formalism is 
compatible with currently employed and established data analysis techniques, and can in principle 
be used to reinterpret existing data in a more general context. The resulting density and pressure 
profiles from this new approach can be used to constrain the equation of state of the galactic fluid, 
and therefore might shed new light on the persistent question of the nature of dark matter. 
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